

' (NASA-CS-1 43862) THE NASTRAN SEB SLAPDOWN N75-25218) 

WATER IMPACT ANALYSIS Final Report 
(Universal Analytics, Inc. , Los Angeles) 

150 p HC $5.75 CSCI 20K Onclas 

V . . _ __ G3/39 24751 / 



Universal Analytics, Inc. 


Los Angeles 



FINAL REPORT 


For 

THE NASTRAN SRB SLAPDOWN 
. WATER IMPACT ANALYSIS 


Presented to: 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
Marshall Space Flight Center, Alabama 


Contract No. NAS8-29665 


May 12, 1975 


Prepared by 


UNIVERSAL ANALYTICS, INC. 
7740 West Manchester Boulevard 
Playa Del Rey, California 90291 

(213) 822-4422 



FOREWORD 


This document presents the Final Report for the contract, "Flexible 
Body/Water Interaction During Water Impact," Contract NAS 8-29 6 65 . The 
project was originated by the Marshall Space Flight Center (MSFC) of the 
National Aeronautics and Space Administration (NASA) , and performed by 
Universal' Analytics, Inc. (UAI) . 

The principal investigator for the project was Mr. D. N. Herting. 
Messrs. R. L. Hoesly, K. Kahyai, E. G. Sergoyan, and D. R. Williamson 
contributed in various phases of the effort. UAI would like to extend 
its appreciation to Mr. D. A. Kross, the NASA contract technical monitor, 
for his able guidance during the performance of the contract and Mr. A. R. 
Thoren of Teledyne Brown Engineering for his assistance in the program 
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ABSTRACT 


This report describes the computer analysis of the hydroelastic inter- 
action between a shell structure representing the Space Shuttle Solid Rocket 
Motor Case and an incompressible fluid during the slapdown phase of water 
impact. The large motions and hydroelastic response of the system is 
obtained by numerical integration of the combined hydrodynamics and structural 
equations of motion. The computerization of the slapdown hydroelastic 
capability has been incorporated into the general purpose NASTRAN Computer 
Code. 

Included in this report are the development of the theoretical basis, 
a guide to the program’s usage, results of correlation and parameter studies, 
and a detailed description of the computer code. 

The effort, performed by Universal Analytics, Inc., under NASA/MSFC 
Contract No. NAS8-29665, is a continuation of the previous task of developing 
water impact loads for the SRB nozzle entry and broadside impact cases. 

Whereas the previous effort treated localized impacts over short time periods, 
the slapdown analysis capability will perform the analysis of large motions 
of the entire booster case for a selected interval after the initial nozzle 
impact. 
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1.0 INTRODUCTION 


The Space Shuttle concept calls for a parachute descent of the expended 
booster solid rocket motors into the sea, followed by recovery. The water 
impact portion of this operation will impart severe loadings which could 
influence the structural design requirements of the solid rocket motors . 
Marshall Space Flight Center (MSFC) has begun a test and analysis program 
to study the dynamic response characteristics of the booster structures and 
define envelopes of acceptable water entry and recovery conditions. 

After the initial impact of the booster nozzle, the vehicle translates 
and rotates into the slapdawn impact condition whereby large pressure loads 
may be applied to 'the cylindrical shell. The slapdown phase of .the space 
shuttle solid rocket booster (SRB) during water entry is a potentially crit- 
ical condition which depends not only on the initial impact condition but 
also on the structural characteristics of the vehicle. . 

This report has been prepared by Universal Analytics, Inc. (UAI) ,to 
describe the theoretical basis, coding design, and results for a computer pro- 
gram to analyze the slapdown impact of a flexible cylindrical structure. The 
slapdown analysis begins after the initial nozzle impact effects have sub- 
sided. The user selects the initial position and velocities at the start of 
the slapdown phase and the total period for the analysis. The computational 
scheme incorporated into NASTRAN analyzes the following effects at 
each instant of time: 

1. The hydrodynamic pressures 

2. The shell response which includes, deformations (velocities and 
accelerations) and stresses in the flexible body 
Vehicle rigid body motion for the total vehicle load 
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The interaction of these functions forms thee coup led hydroelastic analysis. 

The structural model for the impact problem consists of cylindrical finite 
shell elements. The effects of the fluid is modeled by analytically 
calculating the fluid pressures from the structure motions and applying the 
resultant loads directly to the structure. 

The analytic development of the fluid equations is contained in 
Section 2.0. The effects of hydrodynamic pressures, buoyancy, drag, and large 
angle changes in the vehicle- orientation are included. The equations are 
derived with the goal of developing the forces on .the structural finite element 
model which will, in turn, be dependent on the motions of the structural 
degrees of freedom. .The resulting expressions have been cast in a'"form 
conducive to efficient computer processing.- The chapter concludes with a 
discussion of the numerical stability of the transient integration method. 

The method of providing more exact and .stable solutions to the nonlinear itera- 
tions is developed. 

A description of the testing of the computer analysis with finite element 

t 

models of actual structures is given in Section 4.0. Experimental results are 

i 

) 

compared with the analytic values for the 120-inch diameter (77%) test vehicle. 

- * i 

Appendix B contains an extensive set of these comparisons. A full scale SRB 
analytic model was also tested for various impact conditions. ; Also included 
in this chapter are the conclusions and recommendations based on the results 
and the experience with the program. 

The detailed user guide and programming descriptions for the new water 
imapet analysis capability are described in Sections 3.0 and Appendix A 
respectively. The code is implemented in the Level 15.5 NAS TRAN structure 
analysis program as an extension of the existing-'direct transient analysis 
capability. The format of the descriptions, is similar to the format used in 
the NASTRAN User’s and Programmer's Manuals and may be- inserted in the appro- 


priate manuals. 
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2.0 THEORETICAL DEVELOPMENT 


2.1 DEFINITION OF THE SLAPDOWN WATER IMPACT PROBLEM 

When a flexible shell such as the space shuttle solid rocket motor case 
(SRM) impacts upon water, high pressures are generated on the shell surface, 
causing rapid local deformations of the shall and large overall displacements 
and rotations of the center of mass. These deformations and their time rates 
of change interact with the fluid to modify the pattern of flow about the 
shell and, consequently, the pressure profiles and subsequent shell deforma- 
tions. It is essential to accurately predict the forces acting on the vehicle 
and the vehicle’s motion during water impact in order to satisfy the structural 
design requirements of proposed configurations. 

In the slapdown phase of the SRM water impact the following forces act 
on the structure: 

1. The submerged portion of the structure interacts with the x-rater to 
produce a three-dimensional flow field dependent on the vehicle 
position and velocity. The resulting pressures on the shell will 
produce "hydrodynamic” forces which produce local deformations and 
acceleration of the overall rigid body motions. - 

2. Hydrostatic and gravitational forces are produced due to buoyancy 
effects of the submerged volume and mass of the vehicle respectively. 
These forces will primarily contribute to the overall vehicle motion; 
their effect on local deformations is small. 

3. Although the viscous fluid drag forces (i.e., "skin friction") will 

» 

be small when compared to the hydrodynamic forces, they will act as 
damping forces on the overall vehicle motion and affect the results 
when the analysis is performed over the full slapdown time period. 
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These forces are absorbed by the vehicle as a combination of structure 
deformations and inertial reactions, dependent on the structure stiffness and 
mass properties. 

The structure to be analyzed consists of a cylindrical shell with an' 
arbitrary axisymmetric interior structure. The ends of the cylinder are en- 
closed by a nose section and a bulkhead/nozzle structure. The structure may 
be represented by an axisymmetric finite element model whereby the local dis- 
placements at each station are represented by a Fourier series. This allows 
a matrix formulation using uncoupled harmonics and a relatively small number 
of degrees of freedom. Because of the large overall displacements and rota-, 
tions, the dynamics of the structure- will be separated into nonlinear rigid 
body motion and linear local deformations. The deformations will be measured 
in a coordinate system moving with the center of gravity. 

The impact problem must be analyzed by numerical integration of the 
equations of motion with respect to time. At each point in time the. positions , 
velocities, and accelerations of the structure are used to calculate ‘ the cor- 
responding fluid dynamics. The resulting fluid pressures may then be resolved 
into forces on the structure which, in turn, affect the structure motion. 

This numerical transient procedure is developed to achieve accuracy and 
stability within the constraints of economics and computer hardware. Because 
of the large changes in the fluid-structure interface, the fluid forces must 
be treated as nonlinear loads which depend on displacement, velocity, and 
acceleration, rather than as known functions of time. When the nonlinear terms 
become predominant and conventional numerical transient methods are used, a 
numerical stability problem results. In the effort previous to this contract 
[Ref. 1] the standard NASTRAN nonlinear capability proved to have limits on 
the physical parameters of the problem to be solved. 
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In the following sections the theoretical development of. the equations 
of motion and the water impact loads will be presented, followed by a descrip- 
tion of the numerical integration procedure. An overall summary of the 
computational steps concludes the theoretical report. 

1 

2.2 EQUATIONS OF MOTION 

In the following development, the motion of the SRB structure will be 

separated into two types of displacements. The rigid body motion of the 

structure will be measured by the displacements and rotation at the center 

of gravity in a fixed coordinate system shown in Figure 1. In addition, 

the deformations of the structure are assumed to be linear small motions 

measured in a moving coordinate system fixed to the body as shown in Figure 2. 

The structure itself is Assumed to be an axisymmetric shell. At each 

station along the axis, the displacement field is defined by the coefficients 

u n , u?, u 11 of the Fourier series, 
r <p z 
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where u r (<j>) , etc., are the actual displacements at angle, (f>, as shown in 
Figure 2. For a more detailed description of this method of structure model- 
ing, and the options available, see Reference 2. 

In order to separate the local displacements from the rigid body motion 
it is necessary to develop a transformation where:.' 
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FIGURE 1. RIGID BODY COORDINATE SYSTEM 



FIGURE 2. VEHICLE AND ELEMENT COORDINATE SYSTEMS 


( 2 ) 


{ "i } abs ' {u i } local + lD i ] ^ 


CG 


The matrices, D^, in effect, are the free-body modes of the structure. 
Since energy is -preserved in the transformation, the net forces on the center 
of gravity are 


F. 


F > = E [D.] {F.} 


(3) 


M 


1 CG 


When the transformation given above is obtained, the equations of motion may 
be separated into the free body and deflected components as shown in the 
following development. 

The structure is defined by a stiffness matrix [K] , a structural damping 
matrix [B] , and a mass matrix [M] where 

[M] {{i} abs + [B] { " } abs + [K] {u} abs = ' {P} (4) 


where {u} is the vector of absolute displacements at all points on the struc- 
ture and {P} is the corresponding set of applied loads. 

Substituting Eqs. (2) and (3) into Eq. (4), and noting that a rigid body 
displacement does not produce internal loads on the stiffness or damping 
matrices, two separate equations of motion are obtained. 

The equation for local deformation is: 

- M “local + ™ “local + M “local " “ “ { «>C G (5 > 


Since the damping and stiffness matrices do not. restrain rigid body 
motions, the equation for rigid body motion is: 

[D] x [M] [D] {U> cg = [D] T {p} (6) 


Equations (5) and (6) provide the basis for the solution of the problem. 

A set of loads {p} are generated from the fluid pressures, gravity, etc. and 
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transformed to rigid body loads for use in Eq. (6). The rigid body accelerations, 
{U}, are obtained from Eq. (6) and integrated to produce rigid body velocities 
and position. The right-hand side of Eq. (5) is a set of loads on the struc- 
ture which includes the inertial loads due to the rigid body motion. The load 
vector is in equilibrium and, when applied to the structure transient algorithm, 
theoretically will produce no free body motions. 

Note that the matrix [D] is the collection of all partitions [D^] described 
by Eqs. (2) and (3). If the body is undergoing large rotations the [D] matrix 
changes with time, however, the matrix product on the left side of Eq. (6) re- 
mains constant, since the rigid body mass matrix at the CG is defined as: 


[di t [mj m 


M 0 0 
0 M 0 
0 0 I 

— J 


(7) 


where M is the total mass of the structure and I is the inertia. 

Since the matrices [M] , [B] , and [K] are created by the finite element model 
and are developed in Reference 2, Eq. (5) may be solved with the linear struc- 
ture analysis. The remaining tasks are to define the load vector, {P}, and 
the transformation matrix [D] . The latter is derived quite easily in terms 
of the geometry and equilibrium of forces as described below. 

The forces on the structure may be defined as coefficients E n . , F*\ , F n . 

J ri q>x zi 

in the Fourier series representing the .circumferential force distribution at 
each station i.~ The total force -vector in vehicle coordinates at each station 
may be obtained by the equations: 


/ 2lT 

[F ri W cos § ~ sin ^ <3<f> 


F . 
zx 



( 8 ) 
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The total force acting at the center of gravity is the sum of the 
resultant force at each station. The total moment about .the CG is obtained . 
by multiplying the forces by the distance from the center of gravity. 

Rotating the forces into the absolute coordinate system results in the 
following equation for the net force at the center of gravity: 

i 
i 
i 

where 6. is the location of station i relative to the center of gravity. 

i 

The transformation matrices [D^3 defined in Eq. (3) are- the transpose 
of the matrix on the right side of Eq. (9) . In the actual calculation algo- 
rithm these matrices need not be generated explicitly; the loads and acceler- 
ations are transformed as they are needed. 

2.3 VEHICLE LOADS 

In this section the dynamic loads applied to the structure will be 

t 

developed from the position of the vehicle and its displacements, velocities, 
and accelerations . At any point in time the position of the vehicle is 
illustrated in Figure 1. The pressures on each vehicle station are calculated 
individually from the local properties, transformed to loads on the structure, 
and used to calculate the vehicle response. 

2.3.1 Coordinate System Relationships 

At the point of intersection of the vehicle and the water surface at each 
station, two coordinate systems are erected. The types of systems to be used 
in the fluid analysis are shoxim. in Figures 3 and 4. In Figure 3, the basic 
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0 cos 0 
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rigid body positions are used to define the penetration of each section of the 

structure. In Figure 4, a local coordinate system is placed on the surface. to 

define the flow in the region pf each station. 

Referring to Figure 3, the distance from the center line to the water 

surface, r., and the location of the intersection, X., may be obtained from 
a. 3- 

the equations : 

X. = X „ + 6 . cos 9 + r. sin 0 

IL L/Cj 2- 

and . (10) 

r. cos 0 = 6. sin 0 - Z 

i i CCt 

or 

X i = X CG + ^ 6 i " Z CG S±n ^ cos 6 C 11 ) 

2.3.2 Hydroelastic Impact Loads 

The solution for the potential flow problem assuming impact of a cylindri- 
cally shaped body will be summarized in this section. This case is ' represents- 

l 

I 

tive of the broadside slapping condition assuming the booster to be a cylindrical 
shell. Expressions will be developed for three-dimensional pressure ’distri- 

. . I 

butions which will include the interactive effects between the flexible shell 
and the fluid. These loads are calculated for penetration depths which are 
less than the radius of the cylinder. 

Following the procedure described in Reference 1 for determining the 
velocity potential (f>, we shall represent the wetted surface of the penetrating 
shell by an expanding flat plate having a variable width 2a at each vehicle 
station. 

In Figure 4, the coordinate system used for analyzing the potential 
flow is illustrated. The elliptic coordinates n and £ are defined by the 
equations : 
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z 


FIGURE 4. ELLIPSOIDAL COORDINATE SYSTEM AT EACH STATION FOR 

SLAPDOWN IMPACT 
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y = a cosh£ cos n 

. ( 12 ) 

z = a sihhfc, sxn n 

where a is the half-width of the intersection of the vehicle and the fluid 
surface. 

Referring to Figure 4, the half-width, a, of the wetted surface is related 
to the radius of curvature of the shell, R, and to the distance to the water 
line, r_^. If the shell's change in curvature is small in the neighborhood of 
the water line, the following equation results: 

1 

a = (R 2 - r 2 ) 2 , (r < R) (13) . 

i i 

where a and r_^ will vary along the length of the cylinder. 

The use of an elliptic coordinate system is found to be convenient, since 
both the plate- and planar. free surface are distinct coordinate surfaces. The • 
focal distance is equal to the width of the plate. In the elliptic coordinate 
system, Laplace’s Equation is: 


2 - 

V •$ 


(sinh 2 £ + sin 2 n) * 


2 2 
5 § ^ d $ 

2 + 2 


2 

9 $ 

L 2 
3 x 


(14) 


where $ is the velocity potential function. 

A set of solutions to the above equation which conform -to the boundaries- f 
defined in Figure 4 are defined by the equation: 

<p = z(b + B'(X'- X.)"j e _C2m ' i ‘ 1)? sin(2m + l)n' (15) 

m \ m m x / 

where the coefficients B and B', to be determined from the velocities, repre- 

m m 

sent linear variations in the x direction. Note that the choice of a linear 
function for the x direction restricts the velocity parallel to the cylinder 
axis to a constant for each section of the structure. In effect, this is a 
finite difference method for evaluating the effects of flow parallel to the axis. 
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At the moving surface (5 = 0) the coefficients B and B r may be evaluated 


m 


m 


with the following relation between potential and velocity: 




5=0 


ii 

3z 


= - V 


n 


(16) 


and 


3z 


3 $ 35 

35 3z 


1 

a sin f| 35 


(17) 


Equation (15) is evaluated using Eqs. (10) and (17) at 5 = 0. Using the 
orthogonality condition produces the following equations for the coefficients": 


B 

m. 


B’ 

m 


2a 


(2m + 1) 


TT . 

- f v 1 

n .J o n 


i+1 

AB B 1 1 - B 

m _ - m m 

AX 


(q) sin(2m + l)q sinridr) 
i-1 


X. - X. , 
i+l i-l 


(18) 


(19) 


where is the equivalent velocity function at the" surface. - These velocities 
are the sum of both rigid body motion and local deformations. 

Since the known velocities occur 'as Fourier coefficients on the' surface 

! 

l 

of the cylinder it is necessary to develop an approximate method of transfer- 
ring them' to the flat surface, 5 - 0. If V is a velocity normal to the shell 

^ ! 

at an angular position the flow, dQ^, from the shell surface is: ' 


dQ r = V r (<f>) (Rd(f)) (dx cos 6) 


( 20 ) 


The flow through a corresponding section of the flat surface is: 


8Q g - V n dxdy 


( 21 ) 


At this point we introduce a transformation between the points on the shell 
surface and the flat plate. Using a constant one-to-one transformation 


results in the equation: 



Equating the flow through the two surfaces and including the inclination 
effect results in the following velocity transformation: 

V n (n) = f sin" 1 (f) cos e V r c« . (23) 

■The velocity at the shell surface is defined by the Fourier displacement 

coefficients u 11 , where: 
r’ 


N 

v (40 = S cos H't 1 (24) 

r n=0 r 


Substituting Eqs. (23) and (24) into Eq. (18) results in a more 
convenient definition of the flow coefficients : 


where 


and 
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2R sin 1 (“) cos 0 N 
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I 

mn 


*n 
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mn 


• ir 

cos n<f> sin(2m + l)ri sin ri dri 
J 0 " 


4> 


- / a -\ 
sxn ■ (—) 
v r 


cos T) 


(25) 


(26) 


(27) 


Note that the dimensionless integrals are functions only of the width a, 
and may be precomputed numerically as tabular functions. The coefficient 
is obtained from the numerical -difference between two stations. 

The pressure field’ in the fluid may be obtained from the potential 
solution using the linearized Bernoulli equation: 



(28) 


where the derivative of $ must be evaluated at a point fixed in space. The 
value of the potential § at the flat surface- £ = 0 is obtained by substituting 
Eq. (25) into Eq. (16) resulting in: 
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^C-O 


2E I sin 1 (f) 

E E — cos 0 sin (2m + l)ri u 11 

m n . ,v r 

IT (2m + 1) 


(29) 


Since the coordinate system in which we are measuring # is - moving and 
expanding, the following terms are functions of time. 

Because the velocities are measured on the moving body: 


du 11 du 11 

r _ -n . r 
- 7 — = u + — — X 
dt r dx 


(30) 


Because the body is moving relative to the fluid : 


■ &(«^) - 


r. 
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H [6 i SiTl 6 “ Z CG 7 COS 6] 


(31) 

(32) 


(r. tan 0 + <5.)0 — . — 

X x cos 0 


31 


a (obtained by numerical difference) 


mn 


n “ 


_d_ 

dt 


[COS X (J) ]. - f cot 


n 


(33) 

(34) 

(35) 


After performing the appropriate substitutions and different at ions in 
Eqs. (28) and (29), the pressure equation becomes: 


P | = |r PR sin " 1 (J) cos 0 


Z E j [ I (u 11 - i) 

n m ) L' 10 '' r 3x / 


3r 

mn 

. 3a- + / 


^mn \ • - n "1 

TT7) a - J 


n “1 sin ( 2 m + 1)0 

(m + 1) 


+ I (— )u n cos. ( 2 m + l)r] cotr) 
mn a r 


(36) 


The above pressure distributions could be evaluated as a set of points, 
along the flat surface, multiplied by the area to produce forces, and added 
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to the structure. However, since the structure loads must be in the form of 

Fourier coefficients, it will be more convenient to perform a direct trans- 

ttl 

formation from Eq. (36) to the harmonic load, P^, on the K harmonic', where: 



r 


PAX 


2tt 

p((j>) cos k<f> d<j) 


JO 


(k > 0) 


(37) 



r 


RAX 


p(4>) d<j> 


(k -= 0) 


(38) 


In order to conserve energy in the transfer of pressure from the flat 
surface to the shell we observe that if the flow is equal, the pressures 
must be equal, or 


p(<f>) = pI^q 


(39) 


Also, in Eqv. (37), the integration variable, d(f>, may be transformed by 
the equation 


dy = - a sin q dq 


jadcjr 


sxn 


-1 


(“) 

V 


(40) 


dcf> = — sin \ (§) sin q dq 
* K. 


(41) 


* ' * 

- Substituting Eqs . (36) and (41) into Eq. (37) we obtain the k 1 harmonic 

1c 

of the load on station ,i of the surface of the shell, P , where: 

r 


P k = - — ipR 2 AX sin 1 (§) 

r it R 


cos 0 E E 
m n 


r ( -n 3 "r • \ 
\ "r - IF X I 


, _ 31 . -i I I . 

/ 1 mn 1 \ ♦ *n mn mk 

■ ■ mn 3a / R^ - a^ r 2m + 1 



I J . 
mn mk 


(42) 
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and 


ink 


■ TT 

cos k(j> cos (2m + 1) n cos n dp 

. 

0 


(43) 


Note that the same integral functions I ^ are used to transfer loads as were 
used to transfer velocities. The integrals are related complementary 
functions . 

Equation (42) may be used to calculate the fluid impact loads directly. 

The integral functions I and J are precalculated by numerical integration 

for a set of values of a. At each station the values r^, a, a, and x are 

obtained from the rigid body displacements and velocities. If 0 < r . < R and 

a > 0, the velocities and accelerations, u 11 and ii n are obtained from the pre- 

r r 

3 Tm u 

vious solution. The integral values I , J , and — ^ — are obtained by inter- 

& mil’ mn da J 

polating the precalculated values using the local value of a. Equation (42) 
is then evaluated for all harmonics to produce the load vector. 

2.3.3 Additional Loads 

Although the previously developed hydrodynamic impact loads are the 
major cause of structure deformation and internal stresses, the rigid body 
dynamics, required for the above calculations, are also dependent on other 
sources of loading. These loads are the buoyancy, gravity,' and fluid damping 
forces. In addition, the above development is only valid for partially sub- 
merged sections of the structure. Forces due to attached virtual fluid mass 
on the submerged structure will affect the structural response. The above 
effects are described individually in the following discussion. 


Buoyancy Loads 

The buoyancy loads will be added to the '-.dynamic loads at each station. 
At each station, i, the displaced volume is: 
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AV = 


S.., " 6 - 1 

x+1 x-1 


(44) 


where the cross sectional area, S, is: 


S = — 


R_ 

2 


w , .-lx 
TT - 2 sxn — 
K. 


r ± /^ 


Ef - r? 
x 


(45) 


(-R < r. < R) 
x 


•The incremental loads, in structure coordinates, are; 


AP^ = - pgAV cos 8 


(46) 


AP z = PgAV sin 0 


Gravity Loads 

The net gravity load at .the- center of gravity is:. 

P = M • g (47) 

zg 

where M is the total structural mass and g is the acceleration of gravity . 
Viscous Drag Loads 

The viscous drag loads will be applied to the individual structure 
stations in order to provide realistic damping to the structural deformation 
The user provides the coefficients C^ z and The resulting forces at each 

station are: 



where the effective surface area, c., is: 

x 



c . 

X 


(49) 


= r(tt - 2 sin -1 — ) (-R < r ± < R) 

Note that the superscripts 0 and 1 refer to the harmonic and the sub- 
scripts refer to direction of the motion and force. 

Virtual Hass Effects 

When a section of structure is more than one-half submerged, the hydro- 
dynamic impact loads are negligible. However » a certain portion of the 
surrounding fluid remains attached to the vehicle and will modify the struc- 
ture dynamics. The magnitude of the mass may be obtained from the solution 
of two-dimensional flow about a cylinder. A solution to Laplace’s equation 
in cylindrical coordinates is: 

N 

* = E 

n=l 


(at'" + B/ 1 ) 


COS 


n<f) 


(50) 


The boundary conditions are: 


$ = 0 as r °° 

(51) 

V$ « - V = Z u n - cos n<j> at r = E ■ 

r r r 


The potential coefficients 


A 
• n 


and B 

n 


R n+i 


. n. 


are therefore:. 



B = 0 

n 


The pressure function at r = R is therefore: 


P 


= p$ 


P 


R ..n 
— u 
n r 


cos n<f) 


n > 0 


(52) • 


( 53 ) 


The force on each harmonic is defined by the equation: 



p(4>) cos n<|) d(J) n > 0 



r 


RM 

TV 



2ir 


P 


o 

r 


EM 

2tt 



P(*) # 


n = 0 


( 54 ) 


where AH is the length of a section and R is the radius. Substituting 
Eq. (53) into Eq. (54) results in the equation defining the virtual mass on 
each . harmonic n as: 


P 


n 


n 


M u 
nn r 


(56) 


where 


M 

nn 


pR 2 A£ 

n 


n > 0 


57) 


On the zeroth harmonic, a special case, the virtual mass is: 


M = 2pR 2 A& (58) 

oo 

Because these masses change x?hen the vehicle moves into the water, the 
resulting loads, P^, are calculated directly. The masses are scaled accord- 
ing to the wetted area for depths between one-half and fully submerged. 


2.3.4 Loading Application 

The loads P^ and P^ on the structure developed in the previous sections 
are resolved into loads on the rigid body, and the rigid body accelerations 
are. -obtained' (from Eq. ; (6))'.'- Using the- transformations developed -.in Section 2.2 
the inertial reactions are obtained and" added to the structure -load vector, 
resulting in a set of loads that will produce only displacements relative to 
the center of gravity (Eq. (5)). These loads are used in the load vector, N e , 
applied to the structure at each time step. 

The fluid loads are produced by structure accelerations and velocities 
which, in turn, are dependent on the loads. These equations are solved by 
the transient integration algorithm described below which calculates the 
response one time step at a time. 
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2.4 TRANSIENT RESPONSE NUMERICAL METHODS 


In most NASTRAN problems the transient analysis of regular structures 
having predetermined loads is a stable, accurate numerical procedure. However, 
when the loads are functions of the velocities and accelerations , such as the 
water impact of a flexible structure, the numerical stability of the system 
becomes critical. The basic parameters for choosing an integration method are: 

1. Stability : For certain combinations of data, most nonlinear tran- 

sient integration methods may produce a set of errors which grow with 
time until a numerical overflow occurs, indicating a lack of stability 

2. Accuracy : In some cases, the errors remain bounded but eventually 

dominate the solution, causing poor accuracy. 

3. Efficiency : For a given problem size the efficiency of a method 

generally decreases with increased stability and accuracy. Uncon- 
ditionally stable and accurate methods generally tend to require 
very small time steps, many matrix inversions, or' extensive manipu- 
lations for each time step. 

In order to solve the overall problem of hydroelastic fluid/structure 
response during a x^ater impact, a modification of the NASTRAN direct transient 
analysis rigid format has been chosen. 

The basic equation of motion for the water impact analysis is expressed 
by the following matrix equation: 

[M] {u(t) } + [B] {u(t) } + [K] {u(t)} = (P(t)} + (N^(u,u,u) } (59) 

where [M] , [B], and [K] are the structural mass, damping, and stiffness 
matrices, {u(t)} is the vector of displacements of the structural degrees 
of freedom, {P(t)} is a vector of pre-defined loads, and {N„} is the velocity 
and acceleration-dependent fluid load vector which depends on the current state 
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of the system. The method and equations are described in Section 11 of the 
NAS TRAN Theoretical Manual. 


2.4.1 NAS TRAN Algorithm 

In the application of the water impact analysis to the NASTRAN algorithm, 
the loads created by the fluid pressures, including the rigid body forces, 
are entered in the {N^} vector. At each point in time, the displacement, 
velocity, and accelerations of the structure points and CG motions are esti- 
mated from the previous time step. The new loads are calculated using the 
equations given in the previous chapters and an improved set of displacements, 
velocities, and accelerations is computed. The process continues until 
Eq. (59) .is satisfied and the time step is incremented. 

The numerical form of the equations of motion in the NASTRAN algorithm 

is: 


[A 2 ] {u i+1 } 


P. + N„(u.,u.,u.) + [A..] {u.} + A {u. } 

x Z x’ x x l x o x— 1 


(60)' 


where i is the time step index, i.e.: 


{uj = {u(t ± )} 


[a 2 ] = 


At 


1 ,, , 1 . 1 „ 
T M + 2AF B + 3 K 


[A x ] = 


At 


\m-|k 


(61) 


[A q ] = 


At 


1 2 M+ 2S' B -X K 


B i - 3 (u i+l + “i + "i-P 
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and At = t i+ ^ - t^ is the time step size. Note that the solution vector 
{u } is required to calculate the fluid load {M n } since: 

X'T'-L x. 


”i S ^2 <u i+l - 2ll i + “i-P 
4 i " 25F (u i+! - u i-l > 


(62) 


If a vector {u_^ + ^} is found such that Eq. (60) is satisfied it may be 
shown that the HAS TRAN transient integration is unconditionally stable. 
However, if only an approximation from the previous solution for the load, 
is used, the process may diverge as described in Reference 1. 

A method for solving the above problems and correcting for divergent 
characteristics will be developed in the following sections. 


2.4.2 Nonlinear Load Iteration 

At each time step in the transient integration it is desired to obtain 
a displacement vector {u^ + ^} that provides a solution to Eq. (60) . A direct 
method would be to expand the nonlinear loads {N^} in terms of their 
derivatives, i.e.: 


where 


{N^> = - [M f ] {u ± } - [B f ] {u ± > 


M f 

ij 





at u . 
x 


(63) 


(64) 


Substituting the finite difference forms of fu and iL into Eq. (60) we 
obtain: 
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[A 2 ] {u. +1 } - - 


2 M f + 2At E f 
At 


{u i+ i> + {p} 


( 65 ) 


At 


2 M f + A 1 


{u.} + 


1 1 
-±_ M _ ■ x .,. f! 4- A 

2 f 2At f 2 


At 


{u. ..} 

i-l 


Note that the equation has a direct solution if the matrices associated 
with {u.^_^} are combined and inverted. This procedure would be costly, how- 
ever, since a matrix generation and inversion would be required at each time 
step. 

An iterative method may' be used to solve Eq. (65) implicitly whereby 

the solution vector is obtained by repeated substitutions into Eqs. (.62) and 

til 

(60). If y^ is 1 the n approximation to the iterative equations are: 

[a 2 3 {y n+1 > = { N j>,(y n )} + {c> - - . (66) 

where the constant terms are collected in the vector {C} which is: 

{C} = {P ± } + [A-j^] {u.} + [A o ] {u ± _ 1 } (67) 

and the process is started by extrapolating from the last time step by the 
equation: 

{y} = n. + u. , it +^u. . At^ (68) 

o i i-l 2 i-l 

Starting with .the initial estimate, y^, the fluid loads are calculated. ■ 
Equation (66) is solved for the next approximation. The process is repeated 
until the loading error is sufficiently small. The error vector j (e^K-is: 

{e } = {N 0 (y )} + (C) - [A ] {y } (69a) 

n x. n Z n 

or, substituting Eq. (66) with n - 1 for y we obtain: 

. n 

{e „ } " W 1 ’ { Vvi )r ' (69b) 

2.4.3 Stability Analysis 

In some cases the load iteration procedure may be unstable itself. In 
effect i we are solving Eq. (65) in the form: 
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[A - ] {y n+l } ' = ^ [F] {y - } + {C ^ } 


n 


(70) 


where 


m = 


LAt 


■\ M ( + 2S¥ B f 


(71) 


and 

{C 2 > = the constant terms 

Equation (70) is very similar to the equation used in eigenvalue analysis 
whereby the vector contains a set of modes defined by the equation: 


A. [A 2 ] {*.}-- [F] 4.} 


(72) 


The vector, {y }, will be a combination of the exact solution, y , and 
n e 

the modes in the form: 


{y } =-• (y } + £ a? {<}>.} 

n e 3 J 


(73) 


where y^ is the exact solution to Eq. (70), and are the error coefficients. 


Substituting Eq. (73) into Eq. (70) we obtain: 


cu +1 [A 2 ] H.} = - a. [F] {<$>.} (74) 

therefore, comparing Eq. (74) to Eq. (72) we observe that: 

n 1 1 . n / ~r \ 

a. = A. a. ■ (75) 

3 J 1 

If [A [ >1, the error coefficient a_. will grow and the process will 
diverge. Examining Eq. (72) we observe that if the matrix [F] contains larger 
terms than the matrix [A 2 ], some A values will cause divergency. This is very 
possible in a water impact problem when the mass of the attached fluid repre- 
sented' -in the matrix [F] is larger than the structure mass represented in the 

matrix [A„ ’ 

2 
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A similar problem was solved in Reference 3 by the method of extracting 
the actual system modes at each time step using a single step iteration proce- 
dure for a fixed number of modes. This allowed the explicit solution of Eq. 
(59) and supressed the divergent characteristics. In the following discussion 
a similar method will be derived which, in effect, will directly measure and 
remove the errors from the solution. 

2.4.4 Error Elimination Procedure 

In the discussion above it is noted that the load iteration procedure 
produces error vectors that are dominated by the eigenvectors of the matrix 
equation. If a series of the vectors ’ 5 ^n are °^ ta ^ nei ^ by a series of 

iterations, their differences will be rich in the dominant- error modes. As an 
alternative to Eq. (73) , the solution may be expressed as a combination of 
these vectors or: 

{y } = (y } - [U] {a} 

J e. n 

where 

in] - [y x - jy 2 - | ••• j ? n - Vi 1 

The errors are contained in the matrix U, obtained by the iterations, and y^ 
is the exact solution. If the coefficients, a, the weighing factors which, may 
be determined from the equations below, are obtained, the solution is 
available. 

If Eq. (70) is substituted into the basic iteration equation, Eq. (70) , 
the result is : 

[A 2 ] {y n } - [A 2 3 [U] {a} = N(y e ) + {c} (78) 

t 

assuming that the vector {n} is linear for small changes in y, the exact 
load is : 

{N(y e )> = {N(y n )> - [E] [U] {a} (79) 

and from the last iteration: 


:? 6 ) 

(77) 
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therefore 


[F] [U] - [A a ] -tU] 


{a}. '= [A] {y n+1 - y n > 


(81) 


The matrix [F] need not he evaluated explicitly since each column, {U }, 

n 

of [U] may he modified by the identity: * * 


m (u n ) - rn m n - - - [A] {y n+1 - y n } 


(82) 


Combining the above equations results in an expression for the unknown 
coefficients a 


- [A] 


Jl " 


2y i + y o 


! y N+l “ 2y N + y H-l 


{a} = 


EA] {y N+l " V 


(83) 


or 


[A] [A] {a} 


.[A] {6y} 


The above equation results in an exact solution only if the left-hand matrix 
can be inverted. Since the matrix on the left is rectangular it must he trans 

formed to a square matrix. This is performed by pre-multiplying both' sides of 

T 

Eq. (83) by the load matrix U . The coefficients a are approximated by the 
equation: 


{a} 



-1 

[A] [A]l 


[U] T [A] {<5y> 


(84) 


The exact solution is approximated by Eq.' (76) where the above coeffi- 
cients are used directly. This- method will produce good approximations even . 
though the basic iteration equations are diverging. The only requirement for 
satisfactory results is that the number of iterations be larger than the num- 
ber of divergent modes, which, in most cases, is a small number. Also, --note 
that the matrix to be inverted is of the order N-l, where N+l is the number of 
iterations. The matrix is typically small and may be inverted in core. 



3.0 USER’S MANUAL 


3.1 SUMMARY 

This section describes the NASTRAN data deck set up required to execute 
the hydroelastic analysis of a booster vehicle impacting a calm sea. The 
water impact analysis capability in NASTRAN is basically a set of subroutines 
which calculate the fluid loads in a direct transient analysis. The structure 
is modeled with standard NASTRAN elements. The execution of the program 
follows the normal execution path except for the generation of loads on the 
structure. Three impact conditions of the structure are allowed: (a) broad- 

side, (b) tail-first, and (c) slapdown. The finite element idealizations for 
these impact conditions are illustrated in Figures 5 and 6. 

The broadside analysis is used to analyze the impact of a cylinder 
section of structure. This case is a two-dimensional representation of the 
three-dimensional slapdown impact when pressures and local deflections at the 
point of impact are considered. The analysis is valid only for small penetra- 
tion depths . 

The nozzle impact analysis is used for the impact of a conical shell 
structure. Both vertical and limited non-vertical cases are allowed. 
Additional structural elements above the nozzle elements may be included to 
represent the full vehicle dynamics . 

The slapdown option is used to analyze- 'the water impact of a cylindrical 
shell which enters the water at a more general set of initial conditions .' 

The initial height, angle, and velocities may be specified by the user. 
Hydroelastic, hydrostatic, and fluid drag loads are automatically generated. 
The analysis is three-dimensional. Large angular rotation and fluid penetra- 
tion effects are included. 
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SYMMETRY 


Broadside Impact 
(Cylindrical) 



ARBITRARY 

STRUCTURE 


QUAD1 ELEMENTS 


IMPACT 

POINTS 


ORIGINAL PAGE II 
OF POOR QUALIT? 



ARBITRARY 

UPPER 

STRUCTURE 


FIGURE 5. FINITE ELEMENT STRUCTURE MODELS FOR 

BROADSIDE AND TAIL-DOWN CONFIGURATIONS 
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FIGURE 6. FINITE ELEMENT STRUCTURE MODE FOR 
SLAPDOWN CONFIGURATION 


A standard NAS TRAN data deck is constructed with EXECUTIVE CONTROL, 

CASE CONTROL, and BULK DATA sub-declcs, prepared as though a direct transient 
response analysis was requested, with the additions of one case control flag 
(N0NL) and two new bulk data sets (IMPACT and either IMPLIST or SLPLIST) to 
initiate and define the generation of water impact loads. 

3.2 DATA PREPARATION 

This section describes the detailed input data necessary for execution of 
the water impact problems. Features for requesting certain non-standard and 
intermediate output for the impact problems are discussed together with some 
of the more familiar aspects of preparing a NASTRAN transient analysis data 
deck. 

3.2.1 Model Idealization and Executive Control Deck and Preparation 

For the broadside impact case (Figure 5a) the cylindrical booster shell 

\ 

may be modeled in a normal fashion with standard NASTRAN grid points .and BAR 
or QUAD1 elements. The axisymmetric C0NEAX elements must be used exclusively 
for the tail-first nozzle configuration (Figure 5b) or the slapdown case 
(Figure 6) . 

The usual arbitrary selection of grid point locations is allowed, with 
the exception that points which are expected to interact with the fluid 
during the impact process should be spaced uniformly and close together for 
the broadside or nozzle impact models. For the slapdown case, the spacing of 
the elements is arbitrary. 

The modeling tasks of the user include providing the necessary boundary 
conditions and constraints to the structural model. Symmetry in the case of 
broadside impact is defined by modeling one-half of the structure and providing 
single-point constraints on the boundary, points. In the nozzle impact and 
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slapdown cases, the singular rotational displacements must be constrained. 
Samples of input data preparation will be provided in Section 3.3. 

The standard Executive Control data deck is used- with the exception that 
a DIAG = 12 card is interpreted as a flag to delete output of hydrodynamic 
pressures. 


3.2.2 Case Control Deck Preparation 

A typical CASE CONTROL deck includes necessary title information, one 
TSTEP set ID card, and a control flag to exercise the hydroelastic shell/ 
fluid impact analysis (N0NL = 99999). Printed and plotted output of displace- 
ments, velocities, accelerations, non-linear forcing functions, stresses, etc. 
are obtained with standard NASTRAN output request options. Note when axisym- 
metric elements (C0NEAX) are used for the nozzle and slapdown impact analyses, 
grid point labels for output sets are processed by the internal coded format 


ID 


internal 


ID. . „ + 10 
last 


* (n + 1) 


when n = 0, 1, and 2 -for the respective harmonics of response. 


3.2.3 Bulk Data Preparation 

The structural model and transient analysis parameter data defined in 
the Bulk Data deck are input via standard NASTRAN data cards. (See NASTRAN 
User’s Manual.) The properties of the fluid and overall geometric variables 
of the fluid/structure interface are defined with IMPACT and IMPLIST or SLPLIST 
bulk data. The format and data field specification for these cards is given 
in the card description section to follow. In order to provide, an efficient 
and consistent structure/fluid definition, the following special rules are 
imposed for the three impact conditions allowed. 



Cylinder (Broadside) Impact 

a. One-half of the cylindrical shell structure representing the booster 
motor case is modeled. A set of single point- constraints is imposed 
on the vertical axis to represent the symmetric motion. 

b. Planar constraints defined by MPC bulk data cards must be utilized 
when QUAD1 elements are used for modeling. 

c. Physical properties represent a section of the cylindrical shell 
(of length AZ) and dynamic response is assumed localized, at this 
section, i.e., axial interaction effects are ignored. 

d. Evenly spaced grid points (A0 = constant) are required in the region 
of penetration. 

Nozzle (Tail-First) Impact 

a. Axisymmetric (C0NEAX) elements only are allowed for modeling the 
structure, including booster motor case, nozzle, extension, and 
assorted hardware. 

b. Harmonic values 0 and 1, defined on the AX3LC data card, are permitted. 
(0 = vertical entry) . 

c. The canted impact angle (y) between the concial axis and local ver- 
tical must be small and less than the semi-apex angle (a) of’ the cone. 

d. Evenly spaced points (RINGAX) are required along the nozzle in the 
region of penetration. 

Slap down Impact 

a. The structure is modeled with axisymmetric (C0NEAX) elements. The 
portions of the structure not involved with the slapdown loads such 
as the nozzle, nose cone, and internal bulkheads, are also modeled 
with C0NEAX elements. However, these degrees of freedom should be 
partitioned from the matrices via 0MITAX data cards for accuracy and 
more economical processing. 
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b. Up to 20 harmonics are allowed on the AXIC card. The structure dis- 
placements are represented by a Fourier series around the circum- 
ference and the accuracy of the results is dependent on the number of 
terms. However, the running time increases proportionally to the 
number of harmonics. 

c. Although the stations along the structure axis may be arbitrarily 
spaced, it is recommended that wide spacing be avoided. Because 
the pressure is assumed to vary linearly between stations, that 
distance should be less than the radius of the cylinder. 

All Cases 

a. Time step changes in. the integration loop are not permitted. 

b. The iterations performed to refine the water loads within each time 
step are controlled by the convergence criteria parameters ITER, ED, 
and EP. A large value of ITER in conjunction with small ED and EP 
values will improve the accuracy of the results but will increase 
the running time. The use of a smaller time step size with more time 
steps is a less effective alternate. 

3.2.4 IMPACT, IMPLIST, and SLPLIST Bulk Data Descriptions 

As mentioned, new bulk data sets must be supplied when executing the 
hydroelastic water impact analyses. The first set, IMPACT, specifies geo- 
metric and fluid property data uniquely for each of the impact orientations 
of N0ZL, BDSD, or SLAP (i.e., tail-first, broadside, .or - 'slapdown entry, 
respectively). The second set, either IMPLIST or SLPLIST, defines a set of 
grid point identification numbers corresponding to structural grid positions 
on the impacting shell for which hydrodynamic forces may be expected to be 
applied. The SLPLIST card contains, in addition to the grid point numbers, 
their locations along the shell. 
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Definitions of the new bulk data cards are given on the following pages. 
They have been prepared in accordance with the format of the NASTRAN User’s 
Manual and may be directly inserted. 

3.2.5 Slapdown Example Problem 

An example problem is used to illustrate the basic data deck setup and 
program output. The basic problem is described in Figure 7. The NASTRAN data 
deck is shown in Figure 8. Tito C0NEAX elements were used to define a flexible 
cylindrical shell. The three stations on the vehicle are defined by the RINGAX 
data cards. With the exception of the IMPACT and SLPLIST data cards, all of 
the other bulk data are standard NASTRAN cards. The only non-standard NASTRAN 
card in the case control deck is the "N0NL = 99999" card. 

The basic output format is shown in Figures 9 and 10. Figure 9 illustrates 
the printout of the CG motions and the pressure distribution which is output 
for each time step. The position, velocity, and acceleration for each of the 
three rigid body motions is followed by a printout of the pressure at a series 
of circumferential angles' for each submerged station. The station number 
corresponds to the entry number in the SLPLIST data (i.e., the first SLPLIST 
entry is station 1, etc.). Only stations that are less than one-half sub- 
merged produce pressure output. 

Figure 10 illustrates the standard NASTRAN displacement, velocity, or . 
acceleration output as printed versus time. Each particular point corresponds 
to a harmonic number for a RINGFL point as encoded by the equation given in 
Section 3.2.2. Although the element stress output was not requested, it may 
be obtained with the case control request "STRESS =" and with the following 
ALTER cards in the executive control deck: 

ALTER 164,164 

0FP 0EF1,0ES1, , , , // V,N,CARDN0 $ 

ENDALTER 




Water 

(p = .7852 x ICf 4 ) 


' Initial Conditions 
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Cross Section 


FIGURE 7. SLAPDOWN EXAMPLE PROBLEM 


Executive Control Deck 


ID SLAPTEST, EXAMPLE 

APP DISP 

SOL 9,0 

TIME 2 

CEND 

ORIGINAL PAGE IS 
0? POOR QUALITY 


Case Control Deck 


$ 

TITLE = SLAP DOWN EXAMPLE PROBLEM 
LABLE = VV = 100, THETA = 10, HT = 13.0, FLEXIBLE 
$ 

NONLINEAR = 99999 
AXISYM = COSINE 
TSTEP = 20 
SDISP = ALL 
SACCE = ALL 
HARMONICS = ALL 
BEGIN BULK 


Bulk Data Deck 
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SLPLIST 

1 

— 20.0 

2 

-0 

3 20.0 



' TSTEP 

20 

20 
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FIGURE 8 . DATA DECK LISTING FOR EXAMPLE PROBLEM 
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SLAP DOWN EXAMPLE PROBLEM ' * ' FEBRUARY i4V“l975 


'"VV » 100% THETA =“"lO T HT ="l3„0, FLEXIBLE 


TRAN n ENT. S L A P D 6 W N . IMPACT 


o uTp U t ”d A T A 


TIME STEP 14 TIME = 0,06500 

JSENTER OF GRAVI TY P ARAMETERS __ 

HORZ , VERT : ANGLE 
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FIGURE 9 . PRINTOUT FOR CG MOTION AND PRESSURE DISTRIBUTIONS IN SLAPDOWN ANALYSIS 
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SLAP DOWN EXAMPLE PROBLEM 
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FIGURE 10, NAS TRAN PRINTOUT FOR WATER IMPACT ANALYSIS 
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BULK DATA DECE 


Input Data Card IMPACT (Parameters for Water Impact Problem) 

Description : Defines initial conditions, physical constants, and solution 

parameters for transient water impact problems. The basic data 
card has three forms corresponding to the three impact cases . 

FORM 1 (BDSD Case) Broadside Case 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

IMPACT 

TYPE 

RC 

L 

V0 

RH0 

£3 

THMAX 


abc 

IMPACT 

BDSD 

10. 

238. 

240. 

.00012 


20. 


ABC 


+bc 


m 



MM 

ITER 

ED 

EP 


+BC 





8 

fl 

H 

O 

• 

.001 



Field Contents 

RC Radius of cylinder (real, inches) 

XL Length of section of cylinder (real, inches) 

V0 Initial vertical velocity (real, in/sec) 

2 4 

RH0 Mass density of fluid (real, lb-sec /in ) 

THMAX -/Circumferential location of last' point on IMPLIST data (real, 
degrees) 

MM • Number of terms in Fourier series fluid solution (integer, 

default = 10) 


Integration control parameters (see remarks) 

ITER Maximum number of iterations on fluid loads per time step 

(default = 4) 
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ED Convergence criteria on displacement vector, e , magnitude 

difference (real, default: = .01) 

EP Parallel vector test criteria, £p, absolute difference between each 

term (real, default: £p = .001) 

FOBM 2 (N0ZL Case) Inverted Cone Case 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

IMPACT 

TYPE 

m 

RC 

V0 

RH0 

ALPHA 

THMAX 

HMAX | 

abc 

IMPACT 

N0ZL 


100.0 

240. 

.00012 

15.* 

1.0. 

460. : 

1 

ABC 


+bc 

GAMA 

V0L 

PAMB 

0 

NP0LY 

ITER 

ED 

EP 


+BC 

.16 

.1E8 

14.7 


10 

4 

.01 , 

.001 



Field 


Contents • 


RC 

V0 

RH0 

ALPHA 

THMAX 

: -Hl-IAX - 

GAMA 

V0L 

PAMB 

NP0LY 


Radius of base of nozzle (real, inches) 

Initial velocity (real, in/sec) (See Remark #4) 

Mass density of fluid (real, Ib-sec^/in^) 

Half cone angle of nozzle (real, degrees) 

Uniform longitudinal distance between nodes on IMPLIST card 
(real, inches) 

Depth below, water surface for .which, <f> = 0 assumption is made 
(real, inches) 

Impact angle from vertical oblique entry (degrees) 

3 

Initial volume of contained air (real, in ) 

2 

Absolute ambient air pressure (real, lb/in ) 

Number of terras in polynomial fluid solution (integer, default = 5) 


Integration control parameters are the same as for Form 1. 
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FORM 3 (SLAP Case) Slapdown Case 


1 

2 

3 

4 

5 

6 

7 

8 . 

9 

10 

IMPACT 

TYPE 

RC 

THETA 

h ; 

WERT 

VH0R 

VTHETA 

G 

abc 

IMPACT 

SLAP 

10. 

-5. 

i 

j — i 

O 

200. 

0.0 

o 

o 

386.4 

ABC 


+bc 

RH0 

CDR 

CDZ 

PRT 

MM 

ITER 

ED 

EP ' 


+BC 

.001 

0.0 

0.0 

5.0 

4 

4 

.01 

.001 



Field 


Contents 


RC 

THETA 

H 

WERT 

VH0R 

VTHETA 

G 

RH0 

CDR 

CDZ 

NM 

FRT 


Radius of cylinder (real, inches) 

Orientation angle (real, degrees) 

Initial height (from CG to waterline) (real, inches) 

Initial vertical velocity (real, in/sec) 

Initial horizontal velocity (real, in/sec) 

Initial rotational velocity (real, radians /sec) 

Gravitational constant (real, in sec/sec) 

Mass density of fluid (real, lb-sec^/in^) 

Drag coefficient in radial direction (real, dimensionless) 

Drag coefficient in axial direction (real, dimensionless) 

Number of terms in fluid series equations (integer, default = 4) 
Increment at which .pressures are outputed (real, degrees) (only 
if DIAG 12 is not set) 


Integration control parameters are the same as for Form 1, 

Remarks : 

1. The continuation card is required for TYPE = N0ZL and SLAP. The 
default value will be used if a field is blank. 

2. One- and only one IMPACT card is allowed in the bulk data deck. 
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3. 


An IMPLIST bulk data card is required for TYPE = N0ZL or BDSD. 

A SLPLIST (slap list) data card is required for. TYPE = SLAP. 

4. The iteration ITER defines the maximum number of iterations per- 
formed at each time step. The actual number of iterations may be 
less if the actual error is less than the convergence criteria 

e or £ , . When the maximum is reached a closed-form solution is 
P d 

calculated . 

5. If the iteration error on the load vector is less than £ , the 

P 

• vector is assumed to be correct and is used as the final value. 

If the error vectors for two iterations are proportional (i.e., 
only one error mode exists) a closed-form solution is calculated. 

6. The value of the ITER data may be changed to correct for 
instabilities or numerical problems when they occur. It is 
recommended that with small time steps ITER should be decreased. 
With large time steps ITER may be increased for better accuracy. 
Eor any given problem, the optimum time step size is obtained by 
experimentation . 
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BULK DATA DSC* 


Input Data Card IMPLISX (Grid Point Identification - Impact of Shells 
on Fluid) 

Description : Defines the grid point identification numbers for either 

the impact of an inverted cone or horizontal impact of 
cylinder onto a fluid. 


Format and Example: 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

ggig 

SID 

G1 

G2 

• G3 

G4 

G5 

G6 

G7 

+abc 

IMPLISX 

2 

105 

104 

103 

102 

' 101 

100 

99 

ABC 



+abc 

G8 •- . 

G9 

G10 

-etc- 






+BC 

98' 

97 

96 






■ • 

Alternate Form 

* - 


IMPLIST 

SID 

G1 

"THRU" 

GM 






IMPLIST 


105 

THRU 

1 







Field Contents 

SID Grid point set identification number (Integer > 0) 

G1,.,.,GM Grid point identification numbers (Integer > 0) 


Remarks : 

1. In the broadside problem (BDSD) the grid points interacting 

with the fluid must be equally spaced around the circumference, 
starting at the bottom (6=0), with the last' point at location 
0 defined on the IMPACT card. 

"V 


2. 4-13 8 j (1/1/74 - Revised UAI) 
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2 . 


In the iiozzle impact problem (N0ZL) the Gi values refer to 
RINGAX type grid points. The first point lies at the bottom' 
with subsequent points at equally spaced locations along the 
surface. The last point (GM) corresponds to the last node 
which may be submerged in the impact analysis. 

3. The grid-point set ID (SID) is arbitrary and may be used to 
identify a particular run or class of runs. 


2.4-138k (1/1/74 - Revised UAI) 
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BULK DATA DECK 


Input Data Card SLPLIST (Grid Point Identification and Location - Slapdown 
o£ Cylinder on Fluid) 

Description : Defines the grid point identification numbers and axial location 

- for the slapdown of a cylinder on a fluid surface. 

Format and Example : 


'1 23456789 10 


SLPLIST 

Gl 

Z1 

G2 

Z2 

G3 

Z3 

G4 

Z4 

abc 

SLPLIST 

100 

0.0 

1.0 

10. 

102 . 

15. 

103 

25. 

ABC 


+bc 

G5 

Z5 

G6 

Z6 

-etc- 





+BC 

104 

35. 

105 

45. 







Field 


Contents 


G1,...,GM Grid point identification numbers (integer > 0) 
Z1,...,ZM Axial locations of grid points 


Remarks : 

1. The SLPLIST card is used only for the slapdown case. 

2. The grid points must be entered such that the Z values are ascending. 
(Smallest value is Z1 and largest is ZM.) 
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4.0 EXPERIMENTAL CORRELATION AND PARAMETRIC STUDIES 


To measure the validity of the analytic method, accuracy of the program, 
and provide additional insight into the slapdown impact problems, a series 
of computer runs were made to simulate both a test vehicle and the actual 
final design vehicle. For experimental comparison, a substantial amount of 
experimental data was available from drop tests of a 120-inch diameter (77%) 
scale) test vehicle. After correlating with the experimental results, the 
program was applied to a NASTRAN model of the full scale design to determine 
the effects of design differences and parameter changes. This chapter 
summarizes this effort and describes the results. 

Section 4.1 below presents details of the analysis correlation with 
the test. The overall comparison between the test data and the computer 
results was very good and indicates the reasonableness in the engineering ' 
assumptions made in the analysis. 

Parametric studies, performed using both the 77% and full-scale models, 
to test effects of changes in physical’ parameters are presented in 
Section 4.2. This investigation included program executions with various 
values of structural stiffness, structural damping, initial impact position, 
and- initial; velocities .- The- effects - of the. differences between the test 
•vehicle and the- actual structure were also studied. 

Conclusions reached from the experimental correlation and parametric 
studies are given in Section 4.3. 
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4.1 CORRELATION WITH TEST 


The objective of the test program was to provide load and response data 
for derivation of design parameters for the recoverable Solid Rocket Booster 
(SRB) for the Space Shuttle project. The scale model test vehcile, referred 
to as the 120-inch or 77% model, was tested to measure the pressures and loads 
on the vehicle for initial impact, cavity collapse, and slapdown impacts. The 
vehicle consisted of six 120-inch diameter cylindrical sections with a thick- 
ness of 0.375 inches and having the nose, aft bulkhead, and nozzle fabricated 
from existing hardware. 

Because of practical limitations on the instrumentation and occasional 
malfunctions, the experimental data was limited. The instrumentation channels 
had to be divided between the nozzle and the casing. The actual measured 
data from slapdown consisted of 6 diameter deflection gauges, 10 pressure 
transducers, and 54 strain gauges. Further details of the test program are 
given in Reference 4. 

From the analytic viewpoint, the test data was valuable. The measured 
deflections and strains indicated low frequency or nearly static response 
to the changing loads. This provided the practical limits as to time steps 
and frequency content in the analysis. The fluid pressures indicated con- 
centrated loads near the waterline. This indicated that a f ine_ finite - 
— element-' mesh- was-, necessary; - These results guided the’ ■modeling and selection 
of integration parameters for both the scale model and full SRB analyses. 

The NAS TRAN finite element model of the 120-inch diameter test booster 
used for the test comparisons is shown in Figure 11. A listing of the input 
data deck is given in Appendix B. Finite elements, defined with C0NEAX input 
data, were connected to RINGAX points. In this analysis method the displace- 
ments, forces, stresses, etc. are expanded and solved in terms of their 
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FIGURE 11. DEMONSTRATION PROBLEM FINITE ELEMENT MODEL 
120-INCH DIAMETER VEHICLE 





Fourier series coefficients around the circumference. In the structural 
model each harmonic, n, of the Fourier series is uncoupled from the other 
harmonics. This allows economical matrix decomposition and vector solutions. 
Five harmonics were used in the flexible structure analysis (n = 0,1, 2, 3, 4). 
The fifth and higher harmonics were assumed to contribute negligible deforma- 
tion due to their large stiffness. 

The resulting structural mass and stiffness matrices were of order 1410. 
Using structural partitioning, controlled by the 0MITAX data cards, the matri- 
ces were reduced to 741 degrees of freedom. The computer-calculated weight on 
the structure was 90,176 pounds compared to the actual structure weight of 
87,500 pounds. The location of the center of gravity was within two inches - 
of the actual vehicle. In order to restrain spurious structure oscillations 
and simulate the viscous damping of entrapped water, a structural damping 
factor of 5% was added. Although the true damping may be higher, the use of 
the small value was expected to produce conservative results.- 

Initial positions and velocities for the demonstration slapdown problem 

\ 

were taken from the optical trajectory data of drop test C145-020 [Ref. 4]. 

In this test the vehicle was dropped at an initial angle of 30 degrees with 
an initial vertical impact velocity of 52.6 feet per second. Initial condi- 
tions for the slapdown analysis were chosen at a later point in time, when 
the side of the vehicle was just entering the water. The discussion below 
presents typical results from the correlation study. Test data for the 
C145-020 drop was obtained from References 4 and 5; Appendix B contains 
additional correlation data. 

Figure 12 shows the time history of vertical height of the vehicle nose. 
The close comparison of trajectories indicates validity of initial conditions 
and accuracy of resultant loads. 
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Figure 13 presents diameter deflections along the vehicle at 3.0 seconds 
after, initial impact. This time corresponds to maximum deflections at nearly 
all points in the test vehicle; however, the analytical model deflections had 
not yet peaked. The time shift is believed to be due to the bow wave which 
causes the peak pressures to occur sooner. 

Figure 14 shows the envelope of maximum diameter deflections obtained in 
the drop test. The analytical result shows a comparable magnitude (8 inches) 
but occurrence at a location 120 inches forward of the test result. This 
forward displacement is again attributed to the bow wave effect. 

In Figure 15 the pressure time histories at a particular station are 
shown for the centerline (keel), the average over the wetted circumference, 
and the experimental results measured at the keel. The experimental data 
contained considerable noise (+ 10 psi) and a smooth curve was fitted to 
these results. Although the pressures at the centerline were less than the 
test results , the average pressures are generally higher . 

Figure 16 shows pressures versus axial station and circumferential loca- 
tion at 3.1 seconds after impact. Note that with the assumption of a flat 
water surface, the pressures theoretically approach infinity at the waterline. 
In the real case of surface waves, the water is moving away from the vehicle 
at the waterline and the pressures are finite. In the analysis the large 
computed pressures at the waterline compensate for- the smaller computed pres- 
sures at the centerline. Unfortunately, no test data was available for 
pressures off the centerline (keel) . 

The above correlation studies have shown that the analytic results are 
of sufficient quality to use the computer program as a working tool to provide 
analytic backup for the experimental loads program. The overall loads and 
dynamic response and structure deflections appear to match the physical data- 
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FIGURE 15. FLUID PRESSURE COMPARISONS VS. TIME 
120-INCH DIAMETER MODEL 








TIME = 3.1 SEC 




within 20 percent. Primary data differences occurred in the pressure distri- 

i 

butions and the time phase of the loads. Both of these differences may be 
attributed to the effects of surface waves, which cause the following effects: 

1. The existence of a forward bow wave caused the application of pressure 
loads to occur sooner in the test. This was observed in both the 
pressure and deflection results. 

2. The motion of the water along the side of the vehicle probably modi- 
fied the circumferential pressure distributions, causing lower pres- 
sures at the waterline and higher pressures at the keel during the 
actual test. 
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4.2 PARAMETRIC STUDIES 


The objective of the experimental correlation studies was to measure 
the differences between the test' vehicle and the computer results. Analysis 
of the actual SRB structure could then proceed with a known confidence level 
and qualitative understanding of the problem. In conjunction with this 
'analysis, the actual structure will be tested with static loads, obtained 
from experiment for the tested cases . In order to predict the trends of these 
loads and aid in the definition of criteria, additional parametric studies 
were performed on a NASTRAN model of the actual vehicle as discussed 1610X7. 
First, details of the NASTRAN model of the full scale vehicle will be 
presented. 

Illustrated in Figure 17 is the actual rocket booster design, modeled by 
NASA contract support personnel, using the NASTRAN axisymmetric finite ele- 
ments . The fine detail of the nose and nozzle structures x^as included for 
the purpose of analyzing several loading conditions other than the slap down 
water impact. This procedure resulted in a rather, large order analysis with 

I 

3480 degrees of freedom. However, for the slapdoxra. analysis, matrix parti- 
tioning xxras used to reduce the problem to 1033 degrees of freedom. 

The physical .characteristics of the larger booster modified both the 
initial impact of the nozzle and the slapdown dynamics. Although the dia- 
meter was ' increased from 120 inches to . 146- Inches , the overall length went 
.from approximately 1000 inches to a much larger 1535 inches. "This and the 
design of the skirt surrounding the nozzle increased the rotary inertia and 
decreased the angular accelerations. Other significant differences included 
a larger skin thickness of 0.5 inches (vs. 0.375) and stiffer joints betx^een 
the sections. Although the overall size x?as increased, the shell mode 
frequencies also increased due to these increases in stiffness. 
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FIGURE 17. R&OTOTYPE SRB STRUCTURE FINITE ELEMENT MODEL 
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The analysis of the full scale vehicle was performed both by UAI and by 
Teledyne-Brown personnel located near NASA -MSEC. The effort is being contin- 
ued at MSFC. ■ Because of the computer costs involved with the use of a commer- 
cial data center, the. UAI effort was limited. 

The effects of the various physical parameters defining the x-rater impact 
were studied for both the 77% scale model (120-inch diameter) and a computer 
model of the actual SBB (146-inch diameter).. Parameters included shell 
flexibility, damping, and initial conditions. 

Variation of Shell Flexibility 

The 77% structural model was stiffened by truncating the number. of har- 
monics to 0 and 1, which results in no circumferential bending deflection. 

The effect on pressures is shown in Figure 18. In the rigid structure, the 
pressures rise to their peak immediately. In the flexible case, the structure 
motion delays the immediate pressure peak. However, when maximum deflection 
is reached, the pressures and forces may be larger in the flexible case. 

Several analyses of the full SRB model were made using more severe impact 
conditions, scaled from the 77% model tests. An example of the calculated 
pressures for the flexible structure, compared to the rigid case, is shown in 
Figure 19 . The small differences indicate that the pressures on the actual 
vehicle are relatively insensitive to deflections for the stiff full size 
vehicle. 

Effect of Entrapped Water 

Because of the design differences, the full scale vehicle entrapped more 

water than the scale test model. This will cause a substantial change in 

overall weight, as well as increasing damping on the full scale structure. 

The added mass will also lower the angular decelerations and cause higher 

velocities in the forward end of the vehicle. The estimated mass of the 

entrapped water was added to the structure as a fixed distributed mass. 
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FIGURE 18. EFFECTS OF STRUCTURAL STIFFNESS ON PRESSURES 
120-INCH DIAMETER. MODEL 
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FIGURE 19.- EFFECTS ON PRESSURES OF STRUCTURAL STIFFNESS 
146-INCH DIAMETER VEHICLE 
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Variation of Structural Damping 


Aside from the damping effects inherent in the fluid forces, the structure 
and any entrapped water will produce damping. These effects have been observed 
in the oscillations of the experimental data. For the 77% scale model, the 
undamped case is compared in Figure 20 with the results for a structural 
damping ratio of 5%. Although both cases appeared to oscillate more than the 
test data, the damped case more closely resembled the experimental results. 

Effects of Initial Conditions 

In order to investigate the effects of the initial conditions on the full 
SRB vehicle, a series of computer runs were made varying the initial height, 
vertical velocity, and angular velocity at the start of the slapdown process. 
The region of interest is, of course, when the velocities are high and the 
vehicle impacts in a nearly broadside angle. Knowing the effects of the 
different parameters in this region will aid in the future assessment of any 
changes in the vehicle trajectories prior to slapdown. 

Maximum diameter deformation is plotted in Figure 21 versus vertical 
impact velocity for various combinations of height and angular velocities. 

The initial angle was fixed at 27 degrees and the horizontal velocity was 
122 inches per second. Since a change in angle is equivalent to a change in 
height multiplied by the ratio of the angular velocity to the vertical velo- 
city, it was unnecessary to vary the initial angle. 

It was found that horizontal velocity which remains after initial impact 
had little effect upon results. A case having an increase of 50% in horizontal 
velocity was run for this configuration with insignificant changes in diameter 
deformation and CG accelerations. 

Other output quantities such as vehicle accelerations ,. fluid pressures,- 
and structural stresses are also affected by the initial conditions. However, 
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FIGURE 20. EFFECT OF STRUCTURE DAMPING ON FLUID PRESSURES 
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FIGURE 21. EFFECT OF SLAPDOWN INITIAL PARAMETERS ON MAXIMUM DIAMETER DEFLECTIONS 



since the structural response is linear, and the structural vibrations are 
small, the stresses, loads, and thereby the- accelerations are proportional to 
the deformations. For this reason only, the maximum deformation, which indi- 
cates the overall load, is discussed. 

The study of these parameters indicates several interesting effects. The 
first is that the peak loads are not proportional to the square of the initial 
velocities as was expected, but rise more slox^ly. The second is that the 
increase in initial height above the water (or, in effect, a decrease in angle 
of impact) changes the growth rate of the load with velocity. . Both of these 
effects may be explained by the fact that the peak loads occur in the forward 
region of the structure. By this time the impact velocities have been 
decreased by the previous loads on the aft and center of the vehicle. These 
self- cancelling effects insure that the changes in peak stresses and deforma- 
tions will not grow radically with growth in slap down impact velocities. 

In starting the slapdown analysis with the modified NASTRAN program, the 
initial conditions may be selected at any time point in the test trajectory. 

S 

However, because the initial structural deformation is assumed to be .zero, 
it is necessary to start the analysis before any substantial slapdown loads 
and/or deflections have occurred. This point has' generally occurred when 
the angle between the centerline and the water surface is thirty degrees or 
greater. If later initial conditions are specified, structural deflection 
overshoot can occur. Note that this dynamic overshoot only occurs when the 
structural deflections initially have a large value. In some cases, this 
effect was large enough to cause divergence in the numerical procedures 
because the local structural velocities became larger than the overall velo- 
cities. The remedy is to select earlier initial conditions corresponding to 
small deflections. 
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Comparison with 120-Inch Vehicle 

The response of the full scale vehicle was compared to the response of 
the 120-inch test model for the same initial slapdown conditions. If L is 
the ratio of the model length to full scale length, the use of Froude 
scaling results in the following scale factors: 

gravity: G = ~ = 1 

T 

time: T = / L 

velocity: V = — = / L 

1 1 

angular velocity: A = — = 

T /L 

Note that the initial velocity for the full scale vehicle would be larger and 
the angular velocity would be smaller for direct comparisons . 

It was found that the overall scaled pressures and loads in the two 
vehicles were similar, but the full scale model indicated 37% less deflection 
that the scaled test model results. This was due to the fact that the 
increase in skin thickness caused a 135% increase in shell bending stiffness 
(which is proportional to thickness cubed) . This increase in stiffness 
resulted in a smaller coupling effect between the structure deflections and 
the fluid pressures,. 
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4.3 CONCLUSIONS 


The following comments apply to the comparison of the program to the test 
results and to the effects of the various parameters : 

® The correlation studies have shown that peak, overall, xrater loads; 
structure deflections; and vehicle dynamics calculated by the program 
match the test data with error less than 20 percent. 

© The calculated pressures differ from the experimental results both 

in the distribution over the surface and the time at which the initial 
peak occurs. These differences are probably due to the effects of 
surface waves with some additional second order effects of viscousity 
and finite element approximations. 

a The effects of structural flexibility on the resultant fluid loads 
appeared to be small in the 120-inch diameter vehicle and negligible 
for the full scale SRB. 

e The effects of damping due to internal fluid motion and external 

viscous flow are still unknown. These effects are extremely difficult 
to analyze and are not scaled in the small model tests . 
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4'. 4 RECOMMENDATIONS 


The delivered computer program has served its basic purposes of providing 

analytic backup to the experimental testing and measuring the dynamic response 

of the flexible SRB. However, it has proved to be a useful tool for further 

studies of the SRB impact problem, beyond the original intentions and scope 

of the contract. The following are additions to the basic capability which 

could further improve" the capabilities and usefulness . 

o Refine the circumferential pressure distribution by including a 

vorticity correction in the potential function. This would remove 

the theoretically infinite free surface pressure and redistribute 

the circumferential pressure. Keel pressure will increase and 

probably show a better correlation to test values. 

o Provide for longitudinal wave propagation in the potential function. 

This would involve the implementation of additional degrees of 

freedom for the definition of the three-dimensional surface waves: - 

The pressure distribution and the definition of the waterline would 

» 

be modified to account for the pressures and motion at the surface, 
o Miscellaneous improvements could be made to refine the utility of 
the program. These would simplify the user’s tasks of preparing 
the input data and simplify the production of output graphs and . 
summaries. ' These would include additional printout of the waterline 
and diameter deflections at each station, printout of internal error 
criteria for assessing the numerical accuracy, and allowing for the 
definition of input geometry and Initial conditions relative to any 
point on the structure. 
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APPENDIX A - PROGRAMMER 5 S MANUAL 


A.l IMPLEMENTATION INTO NAS TRAN SYSTEM 

The purpose of this appendix is to describe the programming tasks made 
to implement the hydroelastxc water impact analysis into the RASTRAN computer 
program. In this chapter the overall solution algorithm is discussed, in 
A. 2 the modifications to Data Blocks DYNAMICS and NLFT are described, and in 
chapter A. 3 updates to the 1FP, DPD and TED modules are made. The remaining 
chapters describe newly developed fluid data and pressure generating 
subroutines . 

The overall programming scheme has been designed- to provide the user 
maximum convenience within the framework of the existing Level 15.5 NAS TRAN 
System. Modifications and additions to the existing code have been kept to 
a minimum for the following reasons: 

a. Level 15.5 NASTRAN may be updated to a new level and- all changes 
. to existing code will not require reprogramming . 

b. The hydroelastic implementation within NASTRAN will not degrade the 
.system for a normal structural transient analysis. 

c. The reliability and flexibility of the systems is maintained; 
implemented changes involve only isolated subroutines. 

The water impact problem is solved with rigid format 9 (Direct Transient 
Analysis) of the NASTRAN System. The necessary input data given by the user 
are briefly (see User’s Manual for details): 

a. A standard transient analysis structural modal of the booster vehicle 
consisting of axisynunetric shell elements for the nozzle and slapdown 
entries and bar or quadrilateral plate elements for the cylindrical - 
broadside entry. 



b . A list of grid points and coordinate data for the idealized struc- 
ture and points on the- structure entering the water, and 

c. Data defining the fluid properties and various geometric and system 
parameters. 

Except for routine processing of the input data the only differences 
in the program flow from the normal transient analysis lies in the transient 
analysis module (TRD) . The changes to the transient analysis module involve 
replacing 'the operations performed on the non-linear load functions (N0LINI) 
data cards) by the fluid impact operations (IMPACT, IMPLIST, and SLPLIST bulk 
data cards) . 

A simplified flow diagram for the implemented analysis within the TRD 
module is shown in Figure A.l. In order to activate the shell/fluid inter- 
action analysis the control flag N0NLIN = 99999 is set in CASE CONTROL and 
processed in TRD1C. 

A. 2 DATA BLOCK DESCRIPTION UPDATES 

The three new Bulk Data cards IMPACT, IMPLIST, and SLPLIST described in 
the User's Manual (Chapter 3) are processed by’ the IEP module and the card 
images or processed' images written sequentially at the end of the DYNAMICS 
data block. The implementation of these new data- by IEP are described in the 
following pages followed by the DATA BLOCK DESCRIPTION for the input DYNAMICS 
and output NLET tables These pages are prepared In a format suitable for 
direct insertion into the NASTRAN Programmer's Manual 

A. 2.1 Implementation in 'IFF (Table Entries) 

The only code change to subroutine IFP within Level 15.5 NASTRAN consisted 
of a computed G0 T0 'branch to the new bulk data. Implementation required up- 
dates to the IFP block data subroutines IFXiBD (i = 1, 2,..., 7) and subroutine 





FIGURE A. 1 TED MODULE, SIMPLIFIED FLOW DIAGRAM 












IFS3P which was recoded to interrogate the new bulk data cards. These 
changes were : 


IFX1BD: 

IFX2BD: 

IFX3BD: 

IFX4BD: 

IFX5BD: 

IFX6BD: 

IFX7BD: 


Insert IMPACT, IMPLIST, and SLPLIST card ID titles into DATA 16//. 
These occur for card numbers 263 through 265, respectively. 

GIN0 output file 7 and approach acceptability flag 0 set for 
the IMPACT, IMPLIST, and SLPLIST card entries. 

No change 

Required and allowable data items initialized as follows: 

IMPACT 9, 17 
IMPLIST -4, 16 
SLPLIST -4, 9 

(Note: IMPLIST and SLPLIST are open ended cards.) 

Entry statement numbers into IFX7BD field acceptability 
string set as f ollows : 

IMPACT 991 (new string designed to accommodate 16 
non-zero fields) 

IMPLIST -1 (no automated format checking) 

SLPLIST -1 " " " " 

IFP header information consisting of card type identification 
and bit position in a 96-bit trailer assigned as follows: 


IMPACT 

4307, 

43 

IMPLIST 

4107, 

41 

SLPLIST 

4507, 

45 


Format code string entered for unique specification of IMPACT 
bulk data. No automated checking is done on IMPLIST and 
SLPLIST cards. 



DATA BLOCK DESCRIPTIONS 


2. 3. 2. 9 DYNAMICS (TABLE) 


Card Types and Header Information: 



Header Word 

1 

Header Word 2 

Header Word 3 

Card Type 

Card Type 


Trailer Bit Position 

Internal Card Number 

DAREA 

27 

- 

17 

182 

DELAY 

37 


18 

183 

DL0AD 

57 


■ 5 

123 

DPHASE 

77 


19 

184 

EIGB 

107 


1 

86 

EIGC 

207 


2 

87 

EIGP 

257 


4 

158 

EIGR 

307 


3 

85 

EP0INT 

707 


7 

124 

FREQ 

1307 


13 

126 

FREQ! 

1007 


10 

125 

FREQ2 

1107 


11 

166 

IMPACT 

4307 


43 

263 

IMPLIST 

4107 


41 

264 

N0LIN1 

3107 


31 

127 

N0LIN2 

3207 


32 

128 

N0LIN3 ' 

3307 


33 

129 

. N0LIN4 

3407 


34 

-130 

RAMDPS 

2107 


21 

195 

RAN DTI 

2207 


22 

196 

RANDT2* 

2307 


23 

197 

RL0AD1 

'5107 


51 

131 

RL0AD2 

5207 


52 

132 

SEQEP 

5707 


57 

135 

SLPLIST 

4507 


45 

265 

TF 

6207 


62 

136 

TIC 

6607 


66 

137 

TL0AD1 

7107 


71 

138 

TL0AD2 

7207 


72 

139 

TSTEP 

8307 


83 

142 

Type Formats: 

• 


- 


DAREA (4 words) 


SID 

. .• P 

■ C 



A 



DELAY (4 words) 


SID 

P 

c - 



T 



DL0AD (open ended) 

SID 

S 

SI 



L! 

S2 

L2 




S 

L ' 




n 

n 



-1 ' 

-1 


PHASE (4 words) 


SID 

P 

' C 



TH 



EIGB (18 words) 


SID 

METH0D (2 words) 

LI 



L2 

NEP 

NDP 



NON 

E 

N0RM '(2 words) 



G 

C 

0 



0 

0 

0 



0 
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DATA BLOCK AND TABLE DESCRIPTIONS 


Card Type Formats Cont'd.: 


EIGC (open ended) 

SID 

METH0D (2 words) 

N0RM (2 words) 


G 

C 

E 


“a 1 

“a 1 

V 


“bi 


fl e> 


'V 

“a ^ 

V 


V 

V - ' 

. ** 


N e> 

N d 2 

... 


“an 

“an 

“bn 


“bn 

S ' 

•!) 

en 


N dn 


-1 


-1 




-1 . 



EIGP (4 words) 

•SID 

a 



M 



EIGR (18 words) 

SID 

METH0D (2 words) 

-FI 


F2 

NE 

ND 


NZ 

E 

N0RM (2 words) 


G 

C 

0 


0 

0 

0 


0 



EP0INT (1 word) 

ID 



FREQ (open ended) 

SID 

F 

F 


F 

... 

F 


-1 



FREQ1 (4 words) 

SID 

FI 

DF 


NDF 



FREQ2 (4 words) 

SID 

FI 

F2 


... NF 



IMPACT (17 words) 

TYPE (2 words) DATA (12 words, type dependent) 


ITER 

ED 

EP 

IMPLIST (Open Ended) 

SID 

• 6(1) 

6(2) 


G(3) 

6(4) 

6(5) 


♦ , * 

G(M-l) 

G(M) 


-1 



N0LIN1 (8 words) 

SID 

GI 

Cl 


S 

GJ 

CJ 


T 

0 


N0LIN2 (8 words) 

SID 

GI 

Cl 


• S 

GJ 

CJ 


GK 

CK 

- 

N0LIN3 (8 words) 

SID 

GI 

Cl 


S . 

GJ 

CJ 


A 

0 
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Card Type Formats Cont'd.: 


N0LIN4 (8 words) 

SID 

GI 

Cl 


S 

GJ 

CJ 


A 

0 


RAND PS (6 words) 

SID 

0 

K 


X 

Y 

TID 

RANDTI (4 words) 

SID • 

N 

TO 

TMAX 



RANDT2* Not available 




RL0AD1 (6 words) ' • 

SID 

L 

M 


N 

TC 

TD 

RL0AD2 (6 words) 

SID 

L 

' M 

N 

TB 

TP • 

SEQEP (2 words) 

ID 

SEQID 


SLPLIST (open ended) 

G(l) 

z(D 

G(2) 


Z(2) 

etc . 

G(M) 


Z(M) 

-1 

-1 

IF (open ended) 

SID 

GD 

CD 

BO 

B1 

B2 


6(1) 

c(i) 

A0(1) 


Al(l) 

A2(l) 

6(2) 


C(2) 

A0(2) 

AT (2) 


A2(2) 

... 

G(N) 


C(N) 

AO(N) 

A1(N) 


A2(N) 

-1 

-1 


-1 

-1 

-1 

TIC (5 words) 

SID 

G ■ 

C 

1)0 

VO 


TL0AD1 (5 words) 

SID 

L 

M 


0 

TF ‘ 


TL0AD2 (10 words) 

SID 

L 

M 

0 

T1 

T2 


F 

P 

C 


B 



TSTEP (open ended) 

SID 

NO) 

OTO) 

WO) 

N(2) 

DT(2) 


N0(2) 


N(N) 

-1 
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DATA BLOCK AMD TABLE DESCRIPTIONS 


2.3.29.9 FRL (TAB 
Description 

Frequency Response List. 

The FRL contains one logical record for each different frequency set defined in the bulk 
data. Each record contains a sorted list of the unique frequencies defined in the set. 

Table Format 


Record 

Word . 

1X2? 


Item 


•1,2 

B 

Data block 

name 


3 . 

• 

I 

Set ID 

i 

• 



2+n 

I 

♦ 

Set ID„ 
n 



1 

• 

1-m 

R 

Frequencies belonging to set ID ' 

*• 

* 

n 

1-k 

R 

• 

• 

Frequencies belonging to set ID n 

i+l 



End-of-file 


Table Trailer 

' Word 1 ’= number of frequency sets. 

Word 2-6 = zero. 


2.3.29.10 NLFT (TABLE) 


■ Description 

Non-Linear' Forcing Table. 

The header record of the NLFT contains a sorted list of set identification numbers for all 
N0LIN sets defined in the bulk data. Each loqical record of the NLFT contains all data for 
a single set. Point and component numbers on the N0LIN cards are converted to scalar index 
values in both the d- and e-displacement sets. Record 1 in this table is modified 
if a water impact hydrodynamic analysis is desired. 

Table Format 

Record Word Type Item 

0 1,2 B ‘ Data block name 

3 " Set ID 


2+n' Set lD n 
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DATA BLOCK DESCRIPTIONS 


Record Word Type 


Item 


1 


1 

2 

3 

4 

5 

6 


7 


8 


I 

I 

I 

R 

I 

I 


type = 

1 - 

type = 

2 = 

type = 

3 = 

type = 

4 = 

type = 

1 = 

type = 

2 = 

type = 

3 = 

type = 

4 = 


Type of nonlinear load (lstype±4) 

SIL value in d-set 

SIL value in e-set 

Scale factor 

SIL value in d-set 

SIL value in e-set 

Table ID (inteqer) 

SIL value in d-set (inteqer) 
Scale factor (real) 

Scale factor (real) 

Not defined 

SIL value in e-set (inteqer) 

Not defined 
Not defined 


repeated for 
each NHL IN 
card in set 


n 


n+1 


Same format as record 1. 
Data belongs to set ID n . 


End-of-fil e 


Note 


Within each record, the data is sorted on word 2 of each 8-word entry in the record. 


For water impact problems the NLFT output file contains the following 
information: • 


Record 

Word 

Type 

Item 

1 

-L 

2,3, ... ,13 
14, 15, 16 

17 

18 

J J\jXJ 

R and I 
R and I 
I . 

I 

Approach flag BDSD, N0ZL, or SLAP 
Water impact parameters 
Integration control parameters 
Number grid points (ND) 

Number of SIL values (ND0F) 

2 

1,.. .,ND0F 

I 

SIL values in d-set 

3 

1,. . . ,ND 

R 

Z values , axial locations of grid points . 
This record is present only for SLAP case. 

Table Trailer 




Word. 1 = number of N0LIN sets. 


Word 2-6 = zero. 
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A. 3 MODIFICATION TO EXISTING MODULES 


A- 3 -! DPP Module 

In addition to the bulk data processing IFP module described in A. 2.1, 
it is necessary to convert external grid point locations specified on the 
IMPLIST or SLPLIST data card into internal identifiers locating the active 
degree- of- freedom locations in the u d vector set. The following pages describe 
the modifications made to the DPD module for this conversion. These pages are 
prepared in a format suitable for insertion into the NASTRAN Programmer’s 
Manual . 



MODULE FUNCTIONAL DESCRIPTIONS 



f 


o 



1,2,,,., N + 1. 


( 2 ) 


where 



(3} 


3, If the entry corresponds to a FREQ set, the frequencies in the set are sorted, and 
any duplicate frequencies are discarded. The sorted list is written as one logical 
record on the FRL.- 

The RANDPS cards are read into core (if no RANDPS data are present, the PSDL is not 
assembled). The RANDT1 and RANDT2 cards are read into core, and. a list similar to that in 
the frequency processing is formed. This list is sorted on set Identification number. The 
file containing the PSDL is opened to write, and the set identifications are written in the 
header record.' The RANDPS data are written as the first logical record of the PSDL. The 
remainder of the PSDL contains one logical record per set. For RANDT2 sets, the data are 
sorted on time lag, and duplicates are discarded prior to writing the record. For RANDT1 
sets, N + 1 time lags, t-, are written where ' . ‘ * ' 

“ t q + (i-1 )At i 3 1,2 N + 1 (4) 
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MODULE FUNCTIONAL DESCRIPTIONS 


4.47.7.5 Assembly of the NLFT and TRL 


The N0LIN1 (i = 1,2,3, 4) cards are read into core. Each referenced grid point and 
component code is converted to a scalar index value in the u p -set. The data are sorted 
on set identification number. USETD is read into core. The file containing the NLFT data 
block is opened to write, and the set identifications are written in the header record. The 
remainder of the NLFT contains one logical record per set. Scalar index values within each 
set are converted to scalar index values in the and u & sets. The data within each set are 
sorted on the scalar index value to which the forcing function is applied. 


The water impact data cards are processed in a manner similar to the N0LINi 
cards. If the DYNAMICS data block contains records for IMPACT and IMPLIST or 
SLPLIST cards, these data are read into core. The IMPACT data is copies to 
the NLFT file. The list of identification numbers on the IMPLIST (type equals 
BDSD or N0ZL) or SLPLIST (type equals SLAP) are converted to scalar indices in 
the d-set with subroutine DPD4. In general, if the problem is asymmetric 
(TYPE = N0ZL' or SLAP) ; the identification numbers are converted to internal 
numbers by the equation 


ID 


internal 


ID, . . + 10 
list 


* (n + 1) 


where n = 0, 1 for N0ZL and n is user-specified for SLAP case. 

The scalar indices are converted to degree of freedom locations- in the 

• u vector with data block USETD. Each internal point generates two entries: 

u and u~ for the broadside and three entries: u , u , and Ug for the N0ZL and 
r. 0 

SLAP cases. The total number of u^ locations is written on the NLFT file 
followed by the table of u^ pointers (ID0F) . For the SLAP case the list of 
axial locations taken from the SLPLIST card is also written to the NLFT file. 


The TIC cards are read, referenced grid points and component codes are converged to 
scalar index values in the u p -set, and the data are written on SCRl , one logical record per set 
A list of the TIC set identifications is accumulated in core. USETD is read into core. The 
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FUNCTIONAL MODULE DPD (DYNAMICS POOL DISTRIBUTOR) 


file containing the TRL data block Is opened to write, The set identifications are written 
' in the header record. The lasfword of the header contains the degrees of freedom in the u^-set, 
'Data are read from SCSI. Scalar index values are converted to scalar index values in the u^-set. 
Each TIC set is written as one logical record on the TRL. When all the TIC data have been 
processed, the TSTEP data are copied from DYNAMICS to the TRL, one logical record per TSTEP set. 

4,47,7.6 Assembly of the EED and TFL 

Processing of EXG8, EIGC, EISP and E1GR cards is minimal. For each card type present, 
a corresponding logical record is written on EED. For each of the cards which specify P0INT, 
the referenced grid point and component code is converted to a scalar index value (it set ■ 
for EIGB and EIGR cards, u^ set for EIGC cards),' 

Transfer function data -are read from the TF record on DYNAMICS one set at a time. For 
each transfer function set, the point and component codes are converted 'to scalar index values 
in the set, which in turn form row and column numbers of the transfer function matrices. 

The data are written on the TFL, one transfer function set per logical record. The set 
identification number is the first word" of each logical record. Four word entries follow. 

The first word is 65536*coluffin number plus row number* the next three words are the terms of 
the matrices. 

4,47.8 Subroutines 

• Auxiliary subroutines DPD1 , DPD2 , J)pD3, DPD4 are described - above.. 

4,47.8.1 Subroutine Name: DPDAA 

r. Entry Point: DPDAA * , 

- 2, Purpose: To convert a grid point and component; code to a scalar index value In 

the Up set. * 

3, Calling Sequence: CALL DPDAA 

, * 4. Method: ,A flag called INEQ is maintained In /DPDC0M/, If the flag is zero, , - 

EQDYN is read into core and INEQ is set to one. The grid point and component to be . 
converted is stored in BUF(l) and BUF(L-l) where BUF and L are in /DPDC0H/. A binary 
search is performed in EQDYN, If the point is found, the corresponding scalar index 
value is stored in BUF{t), Otherwise, an error message Is queued, and an internal 
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A. 3. 2 TRD Module 


Modifications made to the TRD module for inclusion of the hydroelastic 
analysis occur in several locations and consist of testing for the existence 

of water impact load data (N0NLIN = 99999 in CASE CONTROL) and specific 

calls to the new subroutines. The simplified flow diagram [Figure A.l] 
illustrates the insertion of new subroutines within the TRD module. The 
initialization subroutine, FINTL, reads and converts pertinent bulk data 
(IMPACT, IMPLIST, and SLPLIST) . The routine then sets initial values and 

pointers suitable for one of the three water impact conditions and allocates 

block data and open core for the relative variables and arrays. Subroutine 
NITER performs the solution for each time step using an iterative procedure. 
Maximum number of iterations and convergency criteria are controlled from the 
IMPACT data card. The FPG routine which is called by NITER, is designed 
basically as a driver to either the FC0NE and LGNDR (nozzle entry) , FCYLP 
(broadside cylinder) , or FSLAP and PL0AD (slapdown) subroutines which calculate 
the load vectors at each iteration. The selection of impact type is activated 
by the first word (N0ZL, BDSD, or SLAP) on the IMPACT bulk data card. The 
following pages describe the modifications made to the TRD functional module 
and furnishes descriptions of the new subroutines. They should be inserted 
into the NASTRAN Programmer's Manual, Section 4.65 at the appropriate place. 



FUNCTIONAL MODULE TRD (TRANSIENT ANALYSIS - DISPLACEMENT) 


N0LIN3 loads are computed as. 


P t (t) - 


S {u.(t)}A, u.(t)>0 

J v 


0 , u,(t)<0 

J 


N0LIN4 loads are computed as, 


. -S {~u.(t)}A, u.(t)<0 


P^t) 


0 , u.j(t).>0. 


( 19 ) 


( 20 ) 


For the shell/fluid interaction analysis, the resulting non-linear hydrodynamic 
forces, equivalent to loads processed by the N0LINi bulk data cards, are 
automatically calculated by activating data supplied 'on the IMPACT, IMPLIST, 
and SLPLIST bulk data cards. These non-linear loads are computed by sub- 
routines FC0NE, FCYLP , or FSLAP, having been called through FPG from NITER 
within the time loop in TRD1C. The loads are output directly into open core 
/TRDCl/ for processing the transient analysis solution. The user specifies the 
set of times at which data is to be saved. If the current time is an output 
time, the displacement vector for time t - t^ is output. 

The velocity vector given by: 


b,i - -2St KW - <Vi»> 

is output. 

The acceleration vector given by 

0,1 - [u,.,} - Z {up: 


( 21 ) 


( 22 ) 


is output. 

If the time step is scheduled to change at t^ + -j from At-j to At 2 > the displacement for time 
i+1 has been calculated. {u-}, and (u i+1 > are saved along with (P i+1 l. The matrices 

are formed and decomposed as in Equations 9,10, and 11 for At = 

The following equation is used for computing {u. +2 }> 
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MODULE FUNCTIONAL DESCRIPTIONS 


[D]-{u i+2 } “ j (p] + + P i+2 > + + CC] { u -j + i> + L E J ( u ])* (23) 

The vectors ,{p|) and {u|} in the above equation are calculated as follows. Define 

{u i+l } = ’Kj' ( {u i+l } " {u i } )’ ( 24) 

{ W * :;t ({u i+i } - 2 {u i } + {u i-i })i (25) 

At i 

{G i } = («i+l } " <«W At 2» 

i 

then: ' 

t 2 

[»]) - tu 1+ i> - 4t 2 {i w } 4 (27) 

(28) 

4. Solution of Uncoupled Modal Equation: If the method of matrix formulation is .modal 
and no transfer functions or direct input matrices are used, the equations may be 
solved in a more accurate, more direct manner. The diagonal terms of MUM, BHH, and KUil 
are stored in core. The’ following data are necessary to solve the transient behavior of • 
a modal coordinate (c^)'. 

m^ = Modal mass of mode (MHH) 

.. b.. = Modal damping coefficient (BHH) 

Kj = Modal stiffness (KHH) 



*ot ■ 

(29) 


b 1 

. P i “ 2i7 » 

(30) 


“f ' 3 K\ ~ ^\> 

(31) 


time of the j th time step, 
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A .4 FLUID DATA INITIALIZER (SUBROUTINE FINTL) - 

FUNCTIONAL MODULE TRD (TRANSIENT ANALYSIS - DISPLACEMENT) 

4.65.8.12 Subroutine Name: FINTL 

1. Entry Point: FINTL 

2. Purpose: To initialize shell/fluid interaction data for calculation 

of hydrodynamics loads by the FPG Subroutine. Reads data for conical or 
cylindrical impact problems from NLFT file. Stores pertinent data into 
FPGX common block and initializes pointers fox working storage in open 
core. Allocates open core for use by Subroutine FPG. 

3. Calling Sequence: CALL FINTL (NLFT , IDATA, NDATA, NROW , NGROUP , NZ2 ,IU1 , 

IU2,DELTAT) 

C0MM0N/TRDC1/ - Open core in TRD 

C0MM0N/FP GX/DELZ , DELTH , DT , P I , RC , TL , V0 , RH0 , ALFA ; THMAX , H , GAMA , V0L , PAMB , 

SALFA , CALFA , TALFA , TALFA2 , V0LD , BO , VCR , VC Z , SPY , MM, NP , ND , NK, IT , IPRNT , IPC , 

ICM, ICMS , IP , IRZ , IAS0 , IBS0 , IPMU , IPMUP , IP1M, IP1MU, IAS , IBS , IASD , IBSD , IPZ0 , 
IPZ0P , IP1Z , IP1Z0 , IP0 , IP10 , TYPE , LAST , KHARM, G , CDR , CDZ ,MP , IIP , UD ( 3 ) , UDD (3) , PRT 

NLFT - File number for input data. 

IDATA - Entry point into open core TRDC1 containing water impact 

data and arrays . 

NDATA - Last entry point into open core TRDC1 containing water 

impact data and arrays . 

' NR0W, ) 

NGR0UP f - Pertinent pointers in Subroutine TRD1C. 

NZ2 i 

DELZ - Increment AZ between nodes along longitudinal axis of 

' cylinder or cone. 
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DELTH 


DT 

PI 

RC 

TL 


V0 

RH0 

ALFA 


- Increment A0 between circumferential nodes for cylin- 
drical ' broadside impact, 

- Time increment At in- integration scheme. 

- 3.1415927 

- Radius of cylinder 

- Lentth of cone (for FC0NE) from apex to nozzle exit,’ or 
. length of cylinder (for FCYL 

- Initial impact velocity, V Q . 

- Mass density of fluid 

- Semi-apex angle of cone for nozzle entry 


THMAX 


H 


GAMA 


V0L 


PAMB 


SALFA, 

CA1FA,, 

TALFA, 

TALFA2 


V0LD, 

VCR, VCZ . 


B0 


- Circumferential angle 8 of last point in IMPLIST list 

max r 

for broadside entry (AZ for nozzle entry) . 

- Depth below surface of fluid for which assumption = 0 

» 

holds. Used for nozzle entry only. . 

- Initial canted angle . (y) of nozzle (cone) from vertical 
. (for FC0NE only) . 

- Total volume of enclosed shell impacting vertically in 
inverted conical configuration (for FC0NE only) . ' 

- Ambient air pressure (for FC0NE only) . 

- Sine, cosine, tangent and tangent squared functions of 
semi-apex angle of conical nozzle (for FC0NE only) . 

- Initialized velocity parameters 

- Initialized submerged depth 
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NPY (NP0LY) 
MM 

NP } ND 

NK 

IT 

IPRNT 

IPC - IP 10 

TYPE 

LAST 

KHARM 

G 

CDR 

CDZ 

MP 

IIP 


- Number of terms in expansion of Legendre function for 
stream potential expansion (for EC0NE only). 

- Number of Four'ier terms In expression for stream poten- 
tial (for FCYLP only). 

- Number of nodal points taken along cone or cylinder for 
which hydrodynamic pressures are considered. 

- Number of submerged points on wetted surface for cone or 
cylindrical mathematical models . 

- Time increment counter in integration 'scheme 

- Flag which may be turned on if pressures at each sub- 
merged node and each time step are desired as output 

- Pointers into TKDCl common block for various data and 
working array . 

- Water impact type (N0ZL, BDSD or SLAP) set by first word on 
IMPACT card . 

- Logical value set to true when NITER calls FPG for the 
final time in each time step. 

- Number of harmonics. 

- Gravitational constant. 

- Drag coefficient in radial direction. 

- Drag coefficient in axial direction. 

- Mass property of structure. 

- Inertia property of structure. 
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UD - Rigid body velocities calculated in FSLAP 

TJDD - Rigid body accelerations calculated in FSLAP 

PRT - Increment at which pressures are printed for slapdown 

case (degrees) 

Integration control parameters are passed to NITER in common block IETCRL. 
C0MM0N/INTCRL/ NTER, ED , EP 

NTER - Maximum number of iterations. 

ED - Displacement convergence test. 

EP - Parallel vector convergence test. 

4. Method: The input file (NLFT) is opened and Record 1 is read and 

stored in the appropriate locations in /FPGX/. Work 1 determines the 
approach and the core allocation (N0ZL, BDSD, or SLAP) . See subroutines 
FCYLP, FC0NE, and FSLAP for a description of the core-held data. The 
working data is allocated in open core, words IDATA to NDATA. . The list 
of pointers to other working arrays and vector locations of the’ connected 
grid points (ID0F) is read from the NLFT file and stored in open core. 

For the slapdown case the list of axial locations for each grid point is 
read and stored into open core. 
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A. 5 FLUID PRESSURE GENERATOR SUBROUTINES 

FUNCTIONAL MODULE TED (TRANSIENT ANALYSIS - DISPLACEMENT 

4.65.8.13 Subroutine Name: FPG 

1. Entry Point : FPG 

2. Purpose: To calculate non-linear hydrodynamic loads for shell/fluid 

interaction transient response for either an inverted cone (nozzle) or 
horizontal cylinder impacting a fluid. 

3. Calling Sequence: CALL FPG (DELTAT 5 IC0UNT , IU1,IU2,IP4 > NR0W > LAST,ITP,IY) 

C0MM0N/TRDC1/ - Open core in TRD 

C0MM0N/FPGX/ - See FINTL.. 

DELTAT(DT) - Time increment At of integration scheme. 

IC0UNT(IT) - Iteration counter ' • ' . ' 

IU1 - Index to first location of u.,-, "in TRDC1. 

l+l 

IU2 - Index to first location of u. l0 in TRDG1. 

l+z 

IP 4 - Index to first location of non-linear load vector in TRDCI. 

NR0W - Length of displacement and load vectors. 

LAST - Logical variable to indicate last call to FPG for this 

time step. 

PY - Error flag. Returns integer (1) to FPG if inappropriate 

geometry or loads data is calculated and terminates 
transient response. 

4. Method: FPG is the driver routine which generates the inverted 

concial (FC0NE) , cylindrical (FCYLP) , or slap down (FSLAP) water impact 
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load vectors for transient response. Transmits data pointers into 
TRDC1 for displacements, loads and working arrays. 

4.65.8.14 Subroutine Name: FCYLP (Load Vector Generator for Broadside 

Impact of Cylinder into a Fluid) 

1. Entry Point: FCYLP 

2. Purpose: Calculates non-linear dynamic load vector for the tran- 

sient response of a cylinder impacting a fluid on its side. 

3. '’ Calling Sequence: CALL FCYLP (UD,UDP1,FN,ID0F,CM,CMS,P J IY) 

C0MM0N / TRDC1/ - Open core in TED 

C0MM0N/FPGX/ - Contains variables which, have been initialized in FINTL . 
subroutine- . _ 

UD- - Array in TRD open core containing {u.} — real - input. 

UBPl - Array in TED open core containing {u^_^} - real - input 

■FN - Array in TED open core containing hydrodynamic forces • 

' {ik } real - output and accelerations (u^) as input. 

ID0F - Array of integer pointers. Each entry is the position 

* into the UD, UDP1 and FN vector arrays. 

CM, CMS - Arrays for C at current and past time steps 

.* m 

P - Array of hydrodynamic -pressures 

IY . — Error flag - See FPG 

4. Method 

a. 'Find submerged depth -of penetration, B, and width of wetted 
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B = VAt + B 

previous 


A = / 2R B ' 

c 


where for the first time step B . = 0, V =• V 

previous o. 

b. Obtain velocities at grid points from the structural displace- 
ments, U. , and the ID0F’ table: 

lr 

. ? in - - D ir Ct - 40] / At 

% ' ^(t) - »^(t - AO) / At 

c. Calculate the average, or "rigid body" velocity 



M 



{V. cos 0.'- V., sin 0.} 
in i iq> i 


where M'is the total number of grid points in the IMPLIST list, and 
'0 i = (i - 1)A0. ‘ 

d. Calculate the fluid coefficients C , m = 1, 2...MM using the 
following equations 


*1 





sin 8 . \ 

T- 1 ) 


where ri — ru 
o l l 


IT 


2 and n N+l 


0 
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v. = V. - V cos 9. 
3. m x 


and 


N 

c = /Li V. [sin (2m + 1) uj sin T). T}A. 

m i=l 1 1 11 

where H is the index of the last wetted point . 

e. Calculate the derivative of each C 

. 12 

X ■ .Aft) - C„(t - At)] / At 


The values are stored for use in the next time step, 

f. Calculate the pressures P_^ at each point with the equation: 


P. 

x 


P 

sin u. 
x 


R V • . 2 

c + AV sxn ri . 

A ‘ 


M 



m=o 


1 

2m + 1 



+ A 



sin (2m + 1) n i • 


+ (2m + 1) 

where V = [V(t) - V(t - 
g. The area factors 


R V cos n. 

C — j : — ~~ cos (2m -f 1) rj. 

m A sxn r)- _ x 


At)] /At 


for the points are: 


S. 

x 


- R A9L, 
2 c ’ 


i = 1 
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s. 

1 

R A0L, 
c 

s. 

1 

2 c 


sin 


-1 / A 


(t) 


N-l 


, i = N 


h. The forces' on 'each point in the radial direction are: 

N. = P.S. + u. 

1 1X11 

The forces are placed in the load vector, FN, in positions 
determined by the ID0F table. 


4.65.8.15 Subroutine Name: FC0NE (Load Vector Generator for Inverted Cone 

Impacting into a Fluid) 

1. Entry Point: FC0NE 

2. Purpose: ' To calculate fluid pressure at each submerged node in the 
case of an inverted cone impacting water and to find elements of the 
associated non-linear load vector. 

. 3. Calling Sequence: CALL FC0NE (ID0F , UD , UDP1 , FN , RZ , AS 0 , B S0 , PMU , PMUP , 

P1MU, P1MUP , AS , BS , ASD , BSD ,PZ0 , PZ0P ,P1Z0 , P1Z0P , P0 , PI , IY) 

C0MK0N / TRD C 1 / - Open core in TRD 

C0MM0N/FPGX/ - Contains variables which have been initialized in FINTL 
Subroutine. 


to0F,UD,UDPl 

FN,IY 


See Subroutine FCYLP, each entry of the ID0F array is the 

position -in the UD, UDP1 and FN vectors, e.g. . U n or "u Tl 

n zi 

■for point i and harmonic n. (n = 0 and/or 1) 
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Rz Array of R and 2 coordinates, of nodes on IMPLIST list 

contained in open core /TPDCl/ starting at IRZ. 

AS0 -PI - Working and storage arrays in open core set up in sub- 

routine F-INTL. 

4. Method. Tne logic flow and pertinent mathematical equations are 
as follows: 


a. Penetration depth, B and velocities B and D calculated, and the; 
number of submerged nodes on the cone (NK) determined. 

b. Calculation of geometry (a,^^) of the expanding ellipsoid for 
current time step made. 


•c. Calculate air pressure within entraoped motor case (P ) 

; " a 

d. The velocities (U°) (U°) , (D*) ‘and (D 1 ) ' ' , • . 

r r-t-x z i-rl r x-rl z x+1 are cal-. 


• culated from difference equations for all submerged, nodes' utilizing 

i 

displacements obtained from the UD and UDPi arrays and pointers from 
the ID0F array. ’ ■ ■ 

e. The velocities, normal to the shell, are calculated: 


V n^i^ = CRj) cos « ~ U^Gr) sin 
V n^i^ = ^(ti) cos a - tT^y.) sin 
f- The Legendre polynomials and their first derivative. Namely, 

o *jl o i - 

^m’ ^m an< ^ ^m are calculated for and by calling sub- 
routine LGNDR. 
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g. A and B are calculated as follows, 
s s 


I = 
n 


NK~1 

±=i 


.2 2 




^ IFl o " 

P (]i- ,,) V (U. ,.) -1 
m x+1 n x+1 


^2 P ' 

c - y. 

o x _o , . - , ,, 

—z P (y.) V (y.) 

r 2 _ to x n x 

o 


V i+1 ^1 


TT - 

TO 


NK-l 

Z 

i=l 


\ i p\y ) vj-(y ) + 

s> 


? 2 _ y 2 v 

~P 1 (y .) v^Cy.) 

r 2 ~ i 12 1 n x 

X 


P x+3 y i 


Then: 


A = 
s 


-a (2m + 1) [P;°C? 0 )] 1 1° 


TO 


where m = 2s + 1 


B 

s 


- a 1 ” here ” - 2s + 2 


h. A and B are evaluated as differences with respect to tine 


A - 
s 


A (t) - A (t - At) 

o o 

At 


B 

s 


B (t) - B (t - At) 

b O' 

At 
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i. C , C 91 , C , C 0 „ and C„„ axe calculated for the submerged 


'll* 21’ 12’ 22 

points as follows: 


33 


(1 - p )p • 
_ x i a 

'11 " .2 2 , 

aK o - >V) 


'21 


(C - 1)C 
o o * 

f 2 2. a 

a <L " V 




'12 


Pj. v (i - ft* - l) 
,Jl 2. 


-V sin Y 
o 1 


/ (1 - P.) (C? - 1) 


'22 


a(cj - V J ) 


-V sin y 
o 


C„_ = ~" V sin y 

33 R. o 1 

x 


where V^siny is the average’ lateral velocity of the nozzle. 

■v — v . _ -V* • 

j. P , P , P and P„ are calculated as follows: 

O JL Z o 


NPY 

W - 2- 

s=o u 


A P (U.) P (£ ) + A P ? (p.) P (? ) C 0 , 
s m x m o s m x 21 


+ A P (p.) P‘ (C ) C^ 
s m x m o 11 

J 
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KEY 

VV " 'L, \ v 


s=o 


1 (y.) f 1 ^ ) + a rc 10 p’(y.) ? (& ) 

n i n w o y s [_ 12 iu vl x' m w o' 


+ C P (U ) p»(c >1 

JLJ. to i mol 


+ B 


C-, , p * 1 (y . ) p^(c ) + c 91 p 1 (u.) p ,3 -(c ) i 

11 n x no 21 n x n o 


P 2^i^ 


NPY 

- E \ [c 12 >yV *X> + = 22 p X) lX> 

s=o L -J 


P 3 ^ 


NPY 


- E 


^ B s [ C 3 3 W W] 
s=o L -J 


where 


m = 2s + 1 


n = 2s + 2 


s = 0, 1, 2. . . 


k. The harmonics, n = 0 and 1, of water pressure are e\ r altiated 
for each submerged node from: 


W ” W + f 'W -’'Vh” 


p^Oip - Pj'Cup 


1. Each harmonic of water pressure for an element is assumed to 
■he the average of' that of the adjacent grid points, and is linearly 
distributed to the adjacent grid point in such a manner that the 
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• til 

center of pressure is preserved. The n harmonic of generalized 
normal forces at grid points a and b is: 


an 


F. 


bn 


= 7rL 


= irL 


r a . r b 

% 

L 3 -- 6 -; 

P ) 

av 1 

~r r 1 

a b 

n ( 


P t 

6~ ‘ 3~ 

av ] 


n • ) 


n > 0 


wnere n is the length of the conical element Az/cos a and P„ T is 


av, 


th 


n 


the average of the n harmonic of the water pressure on this • element , 

Note that for n = 0 the forces are multiplied by 2. The forces are 

transformed to the r and z directions, the virtual mass forces M u 

i i 

are calculated, and both are added to the F^ vector in positions 

determined from the ID0F array. The coefficients A and B are saved 

s s 

in locations AS0 and BS0 for use in the next time steo. . 

4.65.8.16 Subroutine Name: LGNDR 

X. Entry Point: LGNDR 

2. Purpose: To calculate the first N odd degree Legendre polynomials 

and first N even degree associated Legendre functions of order one- and 
their first derivatives for variable X. 

3. Calling- Sequence: CALL LGNDR . (N,X,PX,PPX,P1X,PP1X) 

N - Polynomial’s degree (No. of terms in expansion series) 

. integer — input 

X - Given variable - real - input 

PI - Value of polynomial - real - output 
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PPX — Value of derivative of polynomial — real — output 

FIX - Value of function ~ real - output 

PP1X - Value of derivative of function - real - output 

4. Method: The following recursive formulae are used to calculate 

the value of the Legendre polynomials and their derivatives for orders 
0 and 1. 

P*CO - [(2m - 1) XPV^X) - Cm + n ~ 1) P*_ 2 (X) , (tr = 0, 1] 

and 

|~P°(X) = — m [>(X) - P° (X)l 

dxm x 2 _ i j_ m m-1 J 

P 1 (X) = / 1 - X 2 ~ — P° (X) 

m dx m 

and 

4 - P 1 (X) = -ll mXP 1 ^) - (m + 1) P 1 (X) 

dx 2 [_ nr ' v m-1 

x - i 
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Purpose ; To calculate mass and inertia properties at the CG of the structure 
for a slapdown water impact problem. The GG location is also calculated and 
the DLIST array (defined originally on the SLPLIST data card) is converted to 
locations relative to the CG. 

Entry Point : MASSCG 

Calling Sequence ; 

CALL MASSCG (NPTS,NHABM,R,ILIST, DLIST, IM,D,MP,IP,ZC6) 

The variables are: 

NPTS - Number of grid points (input) 

NHARM - Number of harmonics (input) 

R - Radius of cylinder (input) 

ILIST - Array of SIL values (input) 

DLIST - Z locations of grid points (input, output) 
on input - relative to origin 
on output - relative to CG 
IM ~ Matrix control block for MASS file (input) 

D - Working storage (real & integer) 

MP - Mass property (output) 

IP - Inertial property (output) 

ZCG - CG location (output) 

Method : 

The rigid body matrix D is created as a core-held packed matrix. The 
rigid body mass matrix, [M] , is calculated by the equation: 

= [d] t [m] id] 


M 



where M is the structural mass matrix (M^ file) . 


The total mass, m, the CG location, z, and the inertia, X, about the CG 
are calculated by the following equations: 


m = | (M i:l + M 22 ) 


z = - -i M „ 

m 13 


I = - mz 


The location of each station relative to the CG replaces the values z^ 
in the DLIST array, i.e,. 


Z iCG = Z i “ Z 



4.65.8.18 Subroutine Names IJINT, GETIJ 


Purpose : To precalculate the integral terms 1 ^ (a) and. J (a) for a selected 
range of n, m, and a values. The definition of the. integrals is: 

f 1 * 

I -8 cos n<b sin (2m + 1) r| sin p dri 
mxi J q 

r * 

J - l cos nd> cos (2m + 1) p cos r) dp 

mn 

Entry Point : IJINT 

Calling Sequence : 

CALL IJINT (NA,NN,NM,R, Z , ALPHA) 

The variables are: 

Z - Array of open core, output (2*NA*NM*NN) in length 

ALPHA - Array of a values, output 

NA — Number of a points (= -2*MAX(NN/4 NM) + 2) 

NM ~ Number of m values, 1,2,...,NM 

NN - Number of n values, 0,1, — ,NN-1 • 

R - Radius 

Method : 

The routine will generate 2 NA matrices each of dimensions NN by NM and 
store them in the array z. Four nested loops are necessary to perform the 
required calculations. They are: 

1. The outer loop sets the value a. = R0“) and sets the matrices I 

* l NA nm 

and J to zero. The values sin (X. and cos a. are also calculated, 
nm xx 



2. The second loop changes the index k where ri (0 < sin q^ < 1) and 
Aq^ are obtained from the Gaussian integration table. The values 
sin q^, cos n k> sin 2q k , and cos 2q^ are also calculated. 

The functions of <f) are calculated from the equations : 

, cos n 

+ " ““-TTaT 
sin (j) 

The following functions are calculated and stored for 
m - 1,2,. ..,NM. 

CE(m) = cos (2nr + 1) T) = CE(m - 1) cos 2q - SE(m - 1) sin 2q 

SE(m) = sin (2m + 1) q = SE(m - 1) cos 2q - CE(m - 1) sin 2q 

where-. CE(0) =-cos q, and SE(0) = sin q is implied. 

u 1C 

3. The third loop is performed over the index n where n = 0,1,2, ... ,NN. 
The following values are calculated and stored at each value of n > 

CP’(n) = cos (n<j>) = cos (n - 1) <j> cos <f> - sin. (n - 1) -cj> sin eft 

where 

- sin (ncj>) = sin (n -- 1) 4) cos 4> + cos (n - 1) <j> sin cfj 

4. The inner loop is performed over the index m where ra = 1,2,3, ... ,NM. 

The following increments to the integrates are calculated and stored 

Al = Aq, cos nd) sin q sin (2m + 1) q ' 

rat k 

AJ = Aq, cos nrf) cos q cos (2m + 1) q 

nm 'k 


or 



AI - An, CP(n) SE(O) SE(m) 
nm k 

Aj' = An, CP (n) CE(O) CE (m) 

nm k. 

Note that since the indices n and m have zero values, the actual 
program indices are larger by one. 

Entry Point : GETIJ 

Calling Sequence : 

• CALL GETIJ (A,N,M,INM,DINM, JNM) 

The variables are: 

A - Input, real - half width of water-line, a 

N - Input, integer - harmonic index 

M - Input, integer - fluid series index 

INM - Output, real - I integral 

nm 

DINM - Output, real - Sl^/Sa integral 
JNM - Output, real - J - integral 

Method : 

The table is searched until < A < Linear interpolation is 

used for the corresponding values of X and J . The derivative term DINM 

nm nm 

uses the difference AI /Aa. 

nm 



4.65.8.19 Subroutine Name NITER 


Purpose ; To perform the solution of a displacement vector due to nonlinear 
loads for each time step in a transient analysis. 



CALL NITER (DELTAT, IU1, IU2, IU3, IP1, IP4, SCEl, BUF1, IY) 

COMMON /TED Cl/ - Open Core 
The variables are: 

DELTAT - Time step size 

U1,U2,U3 - Pointers to displacement vectors in open core. 

U3 is output 

IPl - Pointer to applied load vector in open core 

IP 4 — Pointer to resultant load vector in open core, output 

SCRljBUEl - GIN0 scratch file and buffer 

IY - Error flag returned by calls to EPG 

Method : 

1. An estimate for the new displacement vector is obtained by linear 
extrapolation, 

{y o > = 2{u 2 > - 

where is the current displacement vector, u^ is the previous 
vector, and y is the estimate of the next vector. 

2. The subroutine performs the iteration on the matrix equation, 

I A 2 ] {y i+1 > = {N( yi )} + {C} 



where: 


[A 2 ] 



1 

2At 



{C} 



(u 2 > + 


M + 
At: 


1 

2At 



{u^ + {P} 


{N( yi )} is the nonlinear load obtained by calling 
Subroutine FPG. 


3- Starting with the second iteration, the following vectors are- 
written to the scratch file for each iteration: 


{6 i } = {y i - y i-i } 


{A i > = {f^ - 


where: 


{f ± } = {N(y ± ) - NCy^)} 


4. The folloxtfing tests are performed after each iteration: 

a. The change in the magnitude of the force vector is checked. 

If jf^l/ |N(y ± ) I < e p , 

the load has converged and the subroutine goes- to Step 6. 

b. The number of non-zero terms in the force vector is checked. 

If i - total number of non-zero terms in {N(y_^)} , 
set n - i and go to 5. 

c. The two previous {f} are checked to ensure they are not parallel.' 

,f . and f. , are normalized such that 
r i-l 




f 


i-1 


= 1 


and the largest term of each is positive. 

If, for all non-zero terms in the vector, 

< e d j = 1,2 NROW 

If all terms pass the test, set n = i - i and go to 5. 

d. The number' of iterations performed is checked. 

If i = I , set n = i and go to 5. 
max ° 

e. If none of the above, continue the iteration for i - i + 1 . 

5. This step is executed if the process is not converging. At least 
three iterations must have occurred. 

a. The equation below is solved for {a}, 

[u] T IA a 3 {a} = [u] T {N (y n ) } 

* 

where the vectors stored on the scratch file make up [u] and 
J A ] as follows : 


h 1 

i — ■ 

ii 

[<$2> 

V 

Ia a ] - 

[’^2 > Ag, . . . 

■ • • • » A^ ] 


Note that the matrix to be inverted is of order ri _ i and may be 
held in core. If the matrix to be inverted is singular, the' 
maximum number of iterations, I max , is decreased and processing 
begins at Step 1. 

b. The final displacement vector is calculated, 

{y} = {y n > “ [«] fa} 


f . . - f . -* . 

f . . 
iJ 



‘c. The final load vector {$[} is calculated by calling 
Subroutine EPG. 

The final displacement and load vectors are placed in the output 
lo.cations indicated by IU3 and IP4, and the subroutine exits. 



4.65.8.20 Subroutine Name FSLAP 


Purpose : To calculate the slapdown loads on the cylindrical shell structure 

for a given set of positions and velocities for slapdown analysis. 


Entry Point: 


FSLAP 



. The common block /FPGX/ contains the control data, miscellaneous con- 
stants, rigid body accelerations and velocities, and total mass properties M 
and I. The displacement vector, U 2 , contains three extra terms corresponding 
to the rigid body displacements U, W, and 0. 

Method : 

The FSLAP subroutine performs three basic steps : The water impact loads 

are calculated for each point for all harmonics and added to the overall 

load vector and to the loads on the CG. The GG accelerations are calculated 

from the basic momentum equations. Finally, the loads on the structure are 

corrected by adding the inertial forces due to the rigid body motion. The 

output load vector, PNL, contains the forces on the structure points plus 

the three total CG forces, F , F , and M n . 

’ x z 0 



1. Estimated displacements and velocities are predicted from initial 
data. The displacements, velocities and accelerations of the 
structure are estimated from the previous time step. 


y i = 


•• . 2 

u. At + 


2u 2i “ u li 


where y is the estimate for displacement, u^ and u £ are the current 
or initial values from the previous time step, and is the esti- 
mated acceleration derived from rigid body considerations. When the 
final iteration is performed, then the last estimate for .accelera- 
tion is calculated as: 


\ = < y i + u li “ 2u 2i )/At ^ 


The detailed processing steps are listed below: 

2. Load Vector Calculations: 

A loop is performed over each station on the structure as 
defined by the table’s ILIST and DLIST. The following steps' are 
executed at each point. 

a. The local velocity functions at point i are: 


.n 


1 / n 

u 

r 


2At (/r - 

.1 If 

1 

1 \ . • 

u r 2At ( 

y r 

- u lr )+» ( 


CG 


i CG (n = 1) 


b. 


where u^ is the displacement at point i in the r direction for 
harmonic n. 

The accelerations are: 


•n 


u 



2u L + 


u lr) 


(n 1) 



u = 


1/1 0 n .1 \ 

77 (y - 2u„ + U, ) 

At 2 V r 2r Ir ) 


+ W CG cos 0 +-U CG sin 0 - 5.0^ 


Note that these equations- are best coded as function subroutines 
to be evaluated as they are needed rather than store all results 
in core. 

c. Calculate the distance from the center line to the water, r^, 
where 

r. = — — — 7 t 5. sin 6 - W p - ui 
l cos 0 [ i CGj 2r 


where u^ is the local displacement at point i in the r direction, 
harmonic 1. 


If > R, the station is "dry", no loads are produced, and 
the loop proceeds to the next point. 

d. If 0 < r < R, perform the following calculations: 


a - 



r . 

i 


(r_^ tan 0 + 6^) 0 


W 
’ * 


cos 0 CG 2At 


(y - u., ) 

r lr 


- a = -• 


r .r . 
x x 


X = U CG + f sin 0 - W CG 6 
Ax = \ (6 i+l “ 6 i-l ) 1 cos 6 


Note 6 = 6' and = 6 >T at the end points 1 and N. 

o 1 N+l N r 

Evaluate the loads P^ using Subroutine PL0AD. The 

61 

integrals I , ~ ~~ and J are obtained by interpolating the 

tables using the current value of a. The derivatives of the 
velocities are obtained from the velocities of the adjoining 


stations : 



5u 


n 


fix 


1 / .n 

2 Ax \ U r 


• n 
— u 


i+1 




where the end values are zero. 

e. The impact loads are distributed to the adjoining stations 
using the equations below: 


F k 

i~l 


_ i P k 

4 r 


E k 

x 


1 pk 

2 r 


1 < i < N. 

x 


F 


i+1 


1 k 
4 r 


The end points on the cylinder are a special case: 


t * 


-|p k 


(i = 1) 


F k = p k =-ip k (i = N) 

N N-l 2 r 

k th 

where is the load in the r direction at point i in the k 

harmonic. The F k forces are added to the appropriate position 

in the PNL vector. 

f. If r, < 0, the loads are calculated from the virtual mass of the 
attached fluid. The net force on each harmonic k is: 


TiT « - “r ' » >0 > 


P° - pMx (R - r; .)) u° (k = 0) 


The loads are distributed to the adjacent stations as in step e. 
If r^ ^ -R, then is set equal to -R-. 



During the execution of step 1 f 1 above, the mass matrix con- 
taining the effective inertia terms relative to the CG is added 
to the running sums M , 6 , and I- . These are added to the mass 
of the structure giving a total rigid body mass matrix, M, where: 



m 
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-M 6 
w w 
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-M 6 
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The new rigid body acceleration in vehicle coordinates is: 



where F , F , and M are the net forces acting at the CG with the 
r z 

acceleration loads removed. 

Inertial Reaction Loads : 

In order to keep the structural displacements small and man- 
ageable, the structure loads are calculated to produce only motions 
relative to the coordinates fixed on the CG. Centrifugal and 
Coriolis forces are ignored. The procedure for calculating these 
corrective loads is : 

a. The mass matrix file is opened to be read a column at a 
time; each column corresponds to a displacement for a given 
point and harmonic. Only the zero and first harmonics are 
used to produce inertial loads. 

b. The loads are generated one point, i, at a time. The. accel- 
erations due to the CG motion at each point are: 
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CG 
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CG 


= u 


"CG 


..1 

u 

z 


= R0 . 


CG 


Note the superscript indicates the harmonic number and the 
subscript indicates the direction, 
c. The columns of the mass matrix corresponding to u^, u^, u°, 
and u^ are multiplied by the negative of the corresponding 
acceleration and added to the load vector PNL in core. This 


load vector will be in equilibrium with respect to rigid 
body motion and will produce only relative structural 
displacements . 



4.65.8.21 Subroutine PL0AD 


j c 

Purpose : To evaluate the loads P for slapdown analysis. 


The equation to b 


be used is: 


,k _ 2 pE . 2 a^ AX cos 6 

r .^2 2 , 

7T (R - a ) 


EE 

n m 


8U 11 

..n r 
u + - 5 ^ 7 - 
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\ a .n 
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u — 
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Entry Point : PL0AD 

Calling Sequence : 

CALL PL0AD (PR,K, MM, XD0T , ADI , UXRD , UKD0T ,DEN , R, AI , DELX,TH,PI, URDD) 
The variables are: 

1c 

PR - Harmonic load (P ) 

r 

K - Number of harmonics • 

XD0T - X 
ADI - a 

.n 

u 

UXRD 

x 

URD0T - u^ 

DEN - p density 

R - Radius 

AX - a 

DELX - AX 

TH - 8 

PI - TT 

URDD - ti 11 
r 

MM - Number of polynomials in series 



Method: 


The routine will precalculate the leading constants. Then two nested 

loops will be used for the summations. Each loop has an entry point for the 

integral subroutine GETIJ. The integral functions I 1 , I , J , are pre- 

me mix tqk. 

calculated and are obtained from tables using the current value of a. The 
following equations are executed inside the loop 


x i - K + 


'2 2 
R - a 


. .n 
X 2 a U r 


r , s, ( ^mn mk 

( *l + x 3 ) 




I j . 

mn mk 


C0W3' = 


2p R Ax cos 9 . -I / a. 


SUM = x. + x c 
4 5 


C0N3 * SUM 



APPENDIX B 


ANALYTICAL AND EXPERIMENTAL DATA COMPARISONS 


Test Run C145-20 


The following plots are given: 

o Height and Angle vs. Time 

© Pressure vs. Time (0=0° and average pressure) 
o Circumferential Pressure Distribution 
o Diameter Deflection vs. Time 

o Diameter Deflection vs. Station 

o Maximum Deflection vs. Station 



Height and Angle vs. Time 
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Pressure vs. Time 


At Centerline and Average Over Wetted Area 


Note: Test -results are plotted as envelopes 

bounded by dashed lines- 
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Run C145-20 

Circumferential Pressure 
Distribution 
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Diameter Deflection 
vs. 

Station 
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Maximum Deflection 
vs . 

Station 







